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Abstract 

Part  I 

Discontinuous  initial  value  problems  for  symmetric  hyperbolic 
lirear  differential  equations 

Initial  value  problems  and  mixed  initial-boundary  value  problems  for 
a  first  order  system  of  partial  differential  equations  are  cctisidered.  Ihe 
initial  data  are  assumed  to  have  jump  discontinuities.  An  existence  and  unique- 
ness theorem  for  weak  solutions  of  such  problems  is  obtained  under  conditions 
which  allow  the  coefficient  matrices  to  have  multiple  eigenvalues,  Th«  proof 
depends  upon  an  a  priori  calculation  of  the  discontinuities  of  the  solution 
and  a  reduction  to  known  theorems  for  problems  with  smooth  initial  data. 

Part  II 

Asymptotic  expansion  of  steady-state  solutions  of  symmetric 
hjrperbolic  linear  differential  eq<iations 

Initial  value  and  mixed  problans  for  the  first  order  system 

Lw  ■  g{x)e     '^     are   considered.     The   solution  is  a  vector  function  w(t,x), 

and  X  denotes  a  vector  of  n  space  variables.     The  vector  g(x)  is  assumed  to 

have  JTiinp  discontinuities.     Under  appropriate   conditions  it  is  shown  that  th« 

solution  has  the  form  w(t,x)  ■  z(t,x)  +  e~      v(x)  where   z(t,x)  -i»  0  as  t  ->co, 

A  formal  procedure  is  introduced  to  obtain  an  expansion  in  powers  of  a>       of 

the  steady-state  function,  v(x) ,     The  results  of  Part  I  are  applied  in  a  proof 

of  the  fact  that  the  fonnal  expansion  is  an  asymptotic  expansion  of  v(x)  for 

CO  ->oo.     The  procedure  and  the  results  are  shown  to  be  applicable  to  Maxwell's 

equations o 
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Part  I 
Discontinuotis  initial  value  problems  for  symmetric  hyperbolic  linear 
differential  equations 

Introduction 

Discontinuous  solutions  of  hyperbolic  partial  differential  equations 
have  often  been  treated  in  the  mathematical  literature.  The  earlier  treat- 
ments dealt  mainly  with  second  order  equations  (see  e.g.  [l] ,  [2]),  Recently 
R.  Courant  and  P.D.  Lax*-  J  extended  the  theory  to  first  order  linear  hyper- 
bolic systems  and  proved  an  existence  theorem  for  discontinuous  solutions 
of  initial  value  problems.  In  Part  I  of  this  paper  we  generalize  the  theory 
of  Courant  and  Lax  in  two  respects.  First  we  consider  more  general  equations. 
This  generalization  is  not  without  interest  for  it  permits  the  inclusion  of 
some  important  equations  of  mathematical  physics.  Second  we  consider  initial- 
boundary  value  problems  rather  than  pure  initial  value  problems. 

In  Part  II  our  theory  will  be  applied  to  the  study  of  asymptotic  expan- 
sions of  solutions  of  certain  partial  differential  equations, 

1,  Initial  value  problems  for  symnetric  hgrperbolic  gelations 

Let  u  be  a  k  component  vector  function  of  the  (n+1)  variables  x  ,  x,, 
,,.,  X  (u  may  be  thou^t  of  as  a  row  or  column  matrix,  according  to  the 
context  in  which  it  appears).  Let  A  and  B  denote  k  x  k  matrices  with  real 
entries.  We  will  sometimes  denote  x  by  t  (since  in  physical  problems  it  may 

V 

represent  time)  and  the  'space  variables'  x,,  ...,  x  by  x.  The  matrices  A 


■? ^ ~ 

That  the  present  theory  is  directly  applicable  to  Maxwell's  equations  is 

demonstrated  in  Part  II, 
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and  B  are  taken  to  be  smooth    functicms  of  x  and  t  in  the  half -space  t  >  0. 
We  define  a  first  order  linear  differential  operator 

(1)  Lu     -     f:    a\     *  Bu,  (u       -  ^  ). 

The  operator  L  will  be  taken  to  be  Bynmetric  hyperbolic,  i.e.,  the  matrices 
A     are  assumed  to  be  symmetric,  and  the  raatrtx  A     to  be  positive -definite. 
These  definitions  are  discussed  in   [U]  and   [5].     Without  loss  of  generality 

we  shall  take  A°  -  I,  where  I  is  the  identity  matrix.     This  is  possible 

o  o    • 

because  the  positive  definite  matrix  A  may  be  written  as  A  ■  T  T,  where 

T  is  a  non-singular  matrix.  Then  the  transformation  of  dependent  vaidables 
Tu  "  V  transforms  L  into  a  synanetric  hyperbolic  operator  on  v  in  which  the 
coefficient  of  v.  is  I . 

Linear  symmetric  hyperbolic  differential  equations  have  been  exten- 
sively studied  in  recent  years  ^  -•  *  ^  -I .  On©  result  of  this  study  is  the 
following  existence,  uniqueness,  and  differentiability  theorem: 

Theorem  It  Let  L  be  a  symmetric  hyperbolic  operator  with  smooth 
coefficients,  0(t,x)  a  smooth  function  defined  in  the  half-space 
t  >  0,  and  g(x)  a  smooth  initial  function  defined  in  the  entire 
x-space.  Then  the  initial  value  problem 

Lu  -  0}  u(0,x)  -  g(x) 

has  a  unique  smooth  solution,  u(t,x), defined  in  the  entire  half- 
space.  The  value  of  the  solution  y  at  a  given  point  depends  on 


Throughout  this  paper,  the  term  'smooth'  will  be  used  to  describe  functions 
which  have  continuous  derivatives  of  sufficiently  high  order  for  the  purpose 
at  hand.  The  precise  order  can  always  be  obtained  by  careful  examination 
of  the  argument,  but  will  usually  be  ignored. 
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the  values  of  the  coefficients,  and  of  the  functions  0  and  g 
only  in  a  finite  region. 

Since  it  is  our  object  to  extend  this  theorem  to  discontinuous  functions, 
we  must  introduce  the  notion  of  a  'weak  solution*.  For  this  purpose  we  con- 
sider a  space  aJ  of  smooth  testing  functions  _Q_(t,x)  with  compact  support* 
More  precisely,  }6  will  be  the  collection  of  all  smooth  k  component  vector 
functions  £\_  each  of  which  vanishes  identically  outside  some  compact  region 
R  which  depends  upon  _Q_  and  lies  in  the  half  space  t  >  0.  For  any  u  in  C  , 
the  set  of  functions  with  continuous  first  derivatives,  and  any  _0  in  jy 
we  have  the  identity 


(2)  _f)_Lu  dxdt  -     u  M  J)_  dxdt. 

Here  M,  the  adjoint  of  L,  is  defined  by 

(3)  Mi\  -  -i:  {/ £\\    -^b'ii. 

V-O  ^       ^  V 

I 

B  denotes  the  transpose  of  the  matrix  B.  The  integration  in  (2)  extends 

over  the  half -space  t  >  0,  The  identity  is  obtained  simply  by  integration 
by  parts,  the  boimdary  terms  vanishing  because  XL  has  compact  support. 
If  now  we  have  Lu  ■  0,  it  follows  trivially  that 

(1a)  J  J  u  M  _fL  dxdt  -  0    for  all  i\  in  s/. 

Conversely  if  u  is  in  C  and  satisfies  (U)  it  follows  that  Lu  »  0  in  the 
half -space,  t  >  0, 

We  may,  however,  consider  more  general  classes  of  functions  than  C  • 

This  region  is  the  'backward  characteristic  cone'.  We  will  not  be  con- 
cerned with  its  precise  delineation. 
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Accordingly  we  introduce  the  class    >(   of  piece-wise  smooth  fxinctionsj 
For  each  vector  function  u  in   A,  defined  in  the  region  of  interest    "^ 
(here   'O    is  the  half -space  t  >  0),  there  are  a  finite  number  of  smooth 
hypersurfaces  which  divide    'C^    into  subregions     /j     such  that  the  first 
derivatives  of  u  can  be  made  continuous  in  the  closure  of  each    X7     by 
appropriate  definition  of  the  derivative  on  the  boundary.       If  a  function 
u  in    A^  satisfies  (U)  we   say  that  u  is  a  weak  solution  of  the  equation 
Lu  ■  0,     If  In  addition  u(0,x)  >  g(x),  we  say  that  u  is  a  weak  solution 
of  the  initial  value  problem 

(5)  Lu  -  0,  u(0,x)  -  g(x). 

We  will  now  analyze  the  discontinuities  and  prove  the  existence  of 
weak  solutions  of  the  problem  (5)  for  appropriate  discontinuous  functions  g, 


2o  The  transport  equations 

Consider  a  function  n  which  is  a  solution  of  (h)  and  which  is  smooth 

everywhere  except  on  a  smooth  hypersurface  C  whose  equation  is  0  (t,x)  -  0, 

u  and  all  of  its  derivatives  are  assumed  to  haw  finite  jump  discontinuities 

across  C»  We  introduce  a  regular  transformation  y.  ■  0''(x^,...,x  ), 

J  •  0,l,«,,,n,  which  utilizes  C  as  a  coordinate  sxirface.  The  equation  of 

C  becomes  y  «  0.  In  the  new  coordinates  we  have 
o 


n 


(6)  Lu 


A 


7*  H*'u   ♦  Bu  • 

3=0      yj 


Here  we  have  introduced  the  matrices  H*'  defined  by 

(7)  H^  -  E  <2'^  A^  •,         j  -  0,1, ...,n. 

V-O    V 


Then  the  discontinuities  of  u  and  its  first  derivatives  will  be  Jump  dis- 
continuities across  the  hypersurfaces. 
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Let  <y  ,  and  /]  ^   denote  regions  on  either  side  of  C  and  let  u  ■  u  +  u 
where  u  i  0  in  -cl  and  u  5  0  in  /o'-.  By  Gauss'  theorem 

(8)         {^-  ^X'^^'^x   d^d*  •  j  j  -^  A^^  ^^^  '  I  (-TLa"^^)  ^If 

w 

Here  a  is  a  scalar  factor  so  chosen  that    (0     ,...,0^   )a  is  a  unit  vector 

o  n 

pointing  in  the  required  direction.     It  follows  that 


(9)  I    I  "^  XLdxdt  -   j   I  XL  Lu'-dxdt  -      f^    (£)_  aV)  a0°  dS. 

C 

2 
By  adding  (9)  to  the  corresponding  equation  for  u  ,  we  obtain  an 

equation  in  which  the  double  integrals  vanish  because  of  (U)  and  the  fact 

that  Lu  ■  0  everywhere  except  on  C.  We  are  left  with 


°  -  .0 


(10)  \Yi    IlA''[u]a0°dS     -  i-LH°[u]adS     -     0. 
\  v^o                                V               )^ 

2  1 

The  synbol  [u]  «  u  -u  denotes  the  Jump  in  u  across  C.  Since  the  components 

of  _f\.  *r«  arbitrary  it  follows  that 

(11)  H°[u]  -  0. 

If  the  vector  [u]  is  not  zero,  the  matrix  H    must  be  singular,  i»e», 

(12)  det  (Yi    0°     A^A  -    0. 

V  v»o      V        / 

But  if  0°  is  a  solution  of  (12),  the  hypersurface  0°  •  0  is  by  definition 
a  characteristic  hypersurface  for  the  operator  L,  Thus  if  u  is  discontinuous 
across  a  hjrpersurf ace  C,  then  C  must  be  characteristic.  We  shall  see  later 
that  the  same  conclusion  follows  if  any  derivative  of  u  is  discontinuous 
across  C. 
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At  this  point  we  inpose  an  important  condition  on  the  class  of 

operators  L  under  consideration,  the 

"  V 

multiplicity  condition:     If  P(t,x,8)  ■  Y_    *yA  ,  where  the  z     are  real 

numbers  not  all  zero,  then  the  multiplicities  of  the  eigenvalues  of  the 

matrix  P  are  independent  of  t,  x,  and  z  ■  (z.,*..,z  ). 

The  multiplicity  condition  assures  us  that  certain  eigenvalues 

(L.,***,a     of  H     are  Identically  zero  on  C  and  the  others  are  non-zero. 

Let  R     denote  the  associated  orthonorraal  eigenvectors.     From  (11),    [u] 

is  in  the  null  space  of  H  ,  and  hence  can  be  expanded  in  terms  of  the 

null  eigenvectors  (i.e.,  those  associated  with  zero  eigenvalues): 

(13)  [u]     -    t    <R"- 

m"l 

The     <r     are  scalars.     Since  Lu  ■  0  on  both  sides  of  0^  -  0, 
(Hi)  E     H^py    1   ♦  b[u]     -     0. 

For  j  f'  0,  u      is  a  derivative  tangential  to  C,  so     u  ■   [u]     •    Thus 

(15)  H^j^yJ   *r    H^Wy     *B[u]     -    0. 

"5 — rp\ ~~~" 

In  the  work  of  Courant  and  LaxL^J,  the  eigenvalues  are  required  to  be  dis- 
tinct (i.e.,   of  multiplicity  one).  For  Maxwell's  equations  P  has  three  eigen- 
values, each  of  multiplicity  two,   (see  Part  II). 

SHI-  O 

It  follows  from  the  lemma  of  Appendix  1  that  the  eigenvalues  of  H  are  con- 
tinuous functions  of  position  on  C.  If  now  cu, ...,a  are  identically  zero 
in  a  region  ^^  on  the  hyperpurface  C,  and  another  set  p, ,,.,,p  were 
identically  zero  in  a  neighboring  region  /(/„  on  G,  then  by  continuity, 
both  sets  would  be  zero  on  the  common  boundary  of  </^  and  /Vp  in 
violation  of  the  multiplicity  condition. 
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Since  H°R^  -  R^°  -  0,  n  »  l,...,q;multipli cation  of  (l5)  by  R^  yields 


(16)  ^   R^'^W       ♦  R^[u]     -0;  [X  '  l,...,q. 


Substituting  (13)  into  (16), 


(17)   r  f  R^^^R"(cr°)  ♦  r  f:  r^^r!  k  *  f:  r^"^™  -  oj  h  - 1 


» 


Equation   (17)  is  a  first  order  linear  syetem  of  partial  differential  equations 
for  the  scalar  coefficients     cT     in  the  expansion  of    [u] ,     It  is  called  the 
transport  equation  of  order  gero«     Before  discussing  it  further,  we  will  de- 
rive the  corresponding  equations  for  the  discontinuities  in  the  derivatives 
of  u.     Front  (6)  we  have  for  r  ■  1,2,*.* 


(18)       ±_  Lu  -  r;     H^u         ♦  Bu       +  r  f"    h£     u   ,     ^v       ♦  N^'^'^^u. 

Here  the  operator  N     Involves  derivatives  of  u  of  order  (r-1)  or  less, 
and  we  have  used  the  abbreviation 

a^-i^u 

u       .   ~- . 

Since  Lu  8  0  on  both  sides  of  C, 

If  all  derivatives  of  u  of  order  less  than  (r+1)  are  continuous  across  C,  then 


A  derivative  of  an  eigenvector  appears  in  (17).  The  lemma  of  Appendix  1 
shows  that  the  eigenvalues  and  eigenvectors  of  H  are  smooth  functions. 
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H°  u  /     ,\|  ■  0.     If  now,  some   (r+l)st,  derivative  has  a  non-sero  jump,  then 

^o 
ju  /     T'vl  /  0  ^^  H    is  singular.     Thvis  we  may  assert  that  if  u  or  any  of  its 

^0 

derivatives  has  a  jiorap  discontinuity  across  C,  then  C  is  a  characteristic 
surface.     Suppose  now  ju       I  is  known  for  s  »  0,l,...,(r-l).     Then  the  jumps 

in  all  derivatives  of  u  of  order  (r-1)  or  less  are  known,  and  hence   JN^^     'u 
is  known.     By  substituting  r  for  (i^l),  we   can  rewrite  (19)  in  the  form 

Here    P^'^     'uj  involves  derivatives  of  u  of  order  (r-1)  or  less,  so  it  is 

known.     We  now  expand     u       I  in  terms  of  all  the  eigenvectors  of  H  J 

^o 

(21)  [u    1  -   r   ,rl  R»  . 

(-y   -i         ffl"l 
"'o 

From  (20)  and  (21)  we  have 

k 


Multiplying  by  r". 

Thus  (23)  gives  the  coefficients  of  the  non-null  eigenvectors  in  the  expansion 

of   lu       I  in  terms  of  known  functions.     To  obtain  the  other  coefficients,  we 

Ly^'J 
^o 

multiply  (19)  by  a  null  eigenvector: 


Any  (r+l)st  derivative  can  be  expressed  as  a  linear  combination  of  u  ^  -,  \  and 
tangential  derivatives  of  r   derivatives  of  u.  Since  all  r       ^o 


derivatives  are  continuous,  their  tangential  derivatives  are  continuous. 


Hence  u  /  ,\  must  be  discontinuous, 
(r+l) 


0 
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l>«*.,q 


o  -J 


It     u  involves  derivatives  of  u  of  order  le?s  than  r,  hence  by  assumption 
is  known.  Upon  substituting  (21)  into  (2I4),  we  obtain 


(25)   f  f  R'-bVc^;)   *  r  f  eVr-^^  ♦  f  R^B.r^  )rV^ 
J"l  ■•!  'J   j"l  m^     'J      m-l       'o 


r 


The  symbol  e    is  defined  by 

(26)   g^;  -  R^^^^^-i^u]  -  r  R^^ 


s«q+l 


<^' 


\i  «  l,...,q 


-  R^BTH°   )  ^     <^^  R% 


^o  s-q+1 


The  g  are  known  functions,  the  (f     being  f^iven  by  (23)  in  terms  of 
lower  order  jumps.  Let  us  establish  the  convention  that  quantities  with 

negative  superscripts  are  zero;  then  since  o^  =  0  for  s  >  q,  g°  -  0  and  (25) 


contains  (17)  as  the  special  case,  r  -  0,  We  rewrite  (25) 


as 


(27)   f  f:  B^R"(0,,*t  ^l.<    ' 
j^  B-l  "^  i       in«l 


g^  ;    i>-  '  i»..«»q 

r  ■  0,1,,,, 


^ere 


(28)     r 


n 


R^^Jr"   ♦  R^(B+rH°  )r"'. 


'i 


It  can  be  sho\«n  that  (2?)  is  for  each  r  a  symmetric  hyperbolic  systen.  There- 
fore we  could  apply  Theorem  1  to  assert  the  existence  and  uniqueness  of  its 

solution.  However  we  can  improve  on  this  approach:  We  shall  show  that  (27) 

* 
reduces  to  a  system  of  ordinary  differential  equations.   To  see  this,  we 

For  this  very  useful  (and  far  from  obvious)  observation,  we  are  indebted 

to  P.  Ungar. 
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shall  prove  the  following 

Lanma:     Let  R  ,.,.,R^  be  orthononnal  eigenvectors  associated 

A-  V     , 

with  the  same  eigenvalue  X    of  the  matrix  P  ■  /"     z  A  ,    (z 

Then  for  each  j. 


(29)  R^  A^  R°     -     c,  6 

0  0  j     tin 

(i«e.,  RqA'^r"     ■     0  if  n  j'  m,  and  is  independent  of  n  if  n  -  ra.) 

Proof:     For  any  j  let  M(K  )  -  P  ♦  ^A^.     Let  R^(i:), .  ..,R^(^)  be  orthononnal 
eigenvectors  of  M  corresponding  to  the  eigenvalue  X(^)  where  X(0)   ■  X     and 
R^(0)  ■  r\     FroB  the  multiplicity  condition,   it  follows  by  «eans  of  the 
lemma  of  Appendix  1  that  the  R  (^)  and  \{k)  are  analytic  functions  in  a 
neighborhood  of  i^  ■  0. 

Let  R^  and  X^  denote  the  derivatives  with  respect  to  J^  evaluated  at 
l^  ■  0  of  R°(^)  and  X(^)  respectively.     Differentiating  the  equation  Mr"  •  Xr" 
we  obtain 

and  multiplying  on  the  left  by  R^, 

R«  A^  r"  -  Ri;(X„-P)R»  ♦  X^rX     •     ^6,™  • 
o  o         ool        jLoo  1  [im 

This  completes  the  proof  of  the  lemma. 

Since  the  matrices  H''  are  just  linear  combinations  of  the  matrices 
A^  it  follows  from  the  lemma  that 

(30)  R^^r"  -  0  If  ii  ^  n',     nyU^R^  -  R^'^R"^}  for  n,ro  -  l,...,q. 
Thus  (27)  nay  be  rewritten 
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<,i,     ^^  rW(^'^)^^  .  t^  tl^  <   .   <  ; 


Expanding  H^  as  in  (7)  yields 


H  ■  l,...,q 

r  »  0,1,.,. 


(32)       f     ^^^^^6^1)       '    t      t    ^i     R^^^R^^^uV. 


V-0  j-1         ^  ^J         V         V"0  ^ 


dy^        n_ 


Here  we  have  used  the  fact  that  ((ff"),     -  0.     Thus   (31)  becomes 


(33)     t   rVrV^),    *t    t^„<    -    g;;;       ^-i,....q 

V"o  V       mil       '^ 

r  ■  0,1,,..     • 

We  now  consider  a  family  of  curves  in  the  tx-space  defined  parametrically 
by  x^  -  x^(8);  V  -  0,1, ...,n}  where 

(3U)  ^v  "  TT  "  ^^^^^^i  v-0,l,..,,n. 

Since  i  ■  £  •  1,  we  may  identify  t  with  s  (up  to  an  ad.iitive  constant  on  each 
curve).  Then  (33)  may  be  interpreted  as  ordinary  differential  equations  along 
these  curves,  for 

t  hVr1( ,.;;),  -r  V^'A  -  K  • 

v«o  V        V»l  "^      V  ^ 


Hence 

(35)  r^  *  t^  rl„<rl  -  g^  j        ^x  -  i,...,q 


nm      n  °|i 

r  .  0,1, 


•  •  •     • 
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We  obsenre  that  ^    ^^  V       '  R""  R  "Oonp^-O,  Hence  the  curves  which 

v«o    *v 
we  have  introduced  are  tangential  to  the  characteristic  surface  0  ■  0»  It 

can  be  shown  that  these  curves,  which  we  shall  refer  to  as  rays,  are  identi- 
cal with  the  characteristic  curves  of  the  single  first  order  partial  differ- 
ential equation  (12)  for  0^,  As  such,  they  are  sometimes  referred  to  as  the 
bicharacteri sties  of  the  operator  L*  From  the  theory  of  first  order  partial 
differential  equations,  we  know  that  the  characteristic  surface  0  •  0  is 
generated  by  an  (n-1)  parameter  family  of  rays* 

Thus  for  each  fixed  r,  (35)  is  a  system  of  q  linear  ordinary  differ- 
ential equations  along  the  rays.  Exairanation  of  the  foregoing  argument 
shows  that  the  coefficients  t'^  and  g'^  will  be  smooth  functions  provided 

y 

the  I'qrpersurface  C  is  smooth  and  the  (T  are  smooth  for  v  <  r.  Given 
appropriate  initial  values  (these  will  be  derived  in  the  next  section),  we 
may  then  assert  that  there  exists  a  unique  smooth  solution  of  (35) •  Thus 
the  system  of  differential  equations  (35)  and  algebraic  equations  (23)  can, 
in  principle,  be  solved  recursively  for  r  ■  0,l,.*.,r.;  and  for  r,  as  large 
as  we  wish.  These  equations  will  be  referred  to  as  the  transport  equations 
(of  order  r). 

Remark:  The  higher  transport  equations  discussed  above  determine  the 
discontinuities  in  derivatives  of  u  normal  to  the  characteristic  surface,  C» 
In  the  special  case  that  C  never  has  a  tangent  plane  parallel  to  the  t-axis, 
we  could  instead  consider  derivatives  in  the  t -direction.  If,  in  addition  , 
the  coefficients  A^,  B  are  independent  of  t  this  leads  to  a  considerable 

simplification  of  the  equations.  A  simplified  analysis  based  on  these 

assumptions  is  carried  out  in  Part  II,  Section  6. 
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3,     Initial  conditions  for  the  transport  equations 

In  the  previous  section  we  were  concerned  with  a  single  characteristic 
surface,  C,  and  its  associated  transport  equations.     Since,  In  what  follows, 
we  shall  have  to  be  concerned  with  several  distinct  characteristic     surfaces 
C^  ,  we  will  use  an  addlticmal  index  (p)  on  various  quantities. 

Let   P  be  a  given  saooth  (n-1)  diaiensianal  manifold  in  the  Initial 
hyperplane,  t  ■  0.     Let  g(x)  be  smooth  on  both  sides  of   p,  while  g  and 
its  derivatives  suffer  juajp  discontinviities  across    p.     Let  Y,,...,Y    be 
components  of  the  unit  normal  to   p ,  and  let  0^     ■  0°^^  (t,x)  -  0  be  the 
equation  of  a  characteristic  surface  C^     whose  intersection  with  the  hyper- 
plane t  -  0  is    p.     Then  0^^^  ■  pT     (v  ■  l,,,.,n)  on   p,  where  p  Is  a  pro- 
portionality factor.*    On    p,  0^^^^   -  ijS^'  is  a  solution  of 

o 

(36)       det/'^    0^P^a'')   -     det  H°^P)  -  det  (0[^^  *  pe)  -    0,  on   p, 
V  v-o     \       / 


idiere 

n 


(37) 


&  ^v  ^- 


Let  R*^i   j  ■  l,...,k  denote  orthonormal  eigenvectors  of  E  (hence  of  H  ^     on 
p)  and  X.  the   corresponding  eigenvalues.     Then  the  eigenvalues  of  H  ^^     ar« 

J 

0+   ♦  P^4«  By  virtue  of  (36),  for  each  p  and  for  some  j, 
t      J 

4  — — — — — .      — _— 

Geometrically,  this  Is  the  condition  that  the  projection  onto  the  hyperplane 

t  «»  0  of  a  vector  normal  to  C^  at  a  point  on  p  shall  be  parallel  to  the 

normal  vector  (Y,,...,Y  )  of  p.  The  condition  is  equivalent  to  the  so-called 

'strip  manifold  condition'  in  the  theory  of  first  order  partial  differential 

equations.  See,  e.g.,  (l].  Chapter  II. 


-lU  . 


(38)  0[P^  +  PX 


There  will  be  as  many  distinct  characteristic  surfaces,  C^  ,  intersecting 
P  as  there  are  distinct  solutions  0|^  of  (38),  hence  as  many  as  there  are 
distinct  eigenvalues  of  E.  The  multiplicity  condition  assures  us  that  the 
nunber  of  these  distinct  characteristic  surfaces  C^*'  is  the  same  at  every 
point  of  p.  On  each  such  C^  a  different  set  of  eigenvectors  is  null  for 
H°^P^  (those  R"^  for  which  X  satisfies  (38)),  but  the  totality  of  null 
eigenvectors  for  all  the  surfaces  C  '^  are  the  totality  of  the  eigenvectors 
of  E. 

Let  Q  ]^,  {_    ]  '  denote  jumps  across  C  '^  and  P  respectively.  Since 
the  eigenvectors  R^  form  a  complete  set  the  Jump  in  initial  values  of  u  across 
P  has  the  unique  decomposition 

(39)  [u]r-  [g]r.f  ^o^R». 

m»l 
We  rewrite  (13),  including  the  additional  index,  p: 


M"  -  r 

m»l 


Here    <r°     «  0  if  r"  is  not  a  null  eigenvector  of  H*'^^   ,     The  initial  conditions 
pm 

for  the    cT       are  obtained  from  the  fact  that 
pm 

(Ul)  E  [u]P  -  M^,     on  p. 

P 

y 

Multiplying  by  any  eigenvector  R  , 
P    ^ 


.  15  - 


But  for  every  v,  R  is  a  null  eigenvector  of  H  on  exactly  one  of  the 


becomes 


characteristic  surfaces,   say  C^''  .     Hence   <r°     -  0  if  p  /  q  and  (U2) 


Thus   (li3)  determines  the  initial  values  of  all  the     cr 

piu 

To  obtain  the  initial  conditions  for  the  higher  transport  equations, 
we  suppose  now  that  lu   \^  have  been  determined  on  P  f°r  ^^  P  ^"^<^  ^°^ 

p  »■  0,1,..,, r-l.  Then  the  jumps  in  all  derivatives  of  order  less  than  r 
can  be  calculated,  (Note  that  we  have  a  different  variable  y  for  each 

of  the  C^^'  and  the  derivative  ~  is  directed  normally  to  C^'^  .)  For 

f    \ 
each  characteristic  surface  C   ,  it  is  convenient  to  introduce  two 

directions  of  differentiation  at  each  point  of  p:  >rr  (normal  to  p  in 

the  plane  t  -  0),  and  -jr-  (normal  to  p  and  tangential  to  C^  ).  The 

corresponding  directions  are  indicated  in  Fig,  1«  Since  the  three 

directions  (y  ,N,s)  are  all  normal  to  P  they  are  linearly  dependent, 

hence  we  may  expand  u   as 

N^ 

r 


(UU)  u   -  H  e^  u  (^_^)  ^ 


N    V"0    y^    s 


Here  the  e^  are  known  angle  factors,  and  e^  /  0  since  the  plane  t  »  0 
is  non -characteristic.   From  (hh) 


A  characteristic  surface  can  never  have  a  tangent  plane  parallel  to  the 


plane  t  «  0.  This  follows  from  (12). 
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Figure    I 


Relative  orientations  of  the  directions  y  ,  N,  s  on  the 
initial  plane. 


-  16  - 


UN  -J      i-y, -"    v«l    '-yl   -Js 


y  -J    v«i    '— y 
J'o  •'o 


By  (U5)>  ju   I  and  u      daterroine  each  other  uniquely.  As  before,  we 

o 


rewrite  (21)  as 

()46)  fu    1^    -   I:    crl^   r" 

Let 

k 

'     L 
.     pro 


(U7)         |u    1^  .  r  pL  r"  . 


[>]' 


Then  from 

(U5) 

(U8) 

where 

(U9) 

<. 

^tvT  ■  •-  "-•  *  ^ 


o         pn  pix 


^^  L  ^v  L(.-v,] 


p 

V 


By  (U8),  the  B       and   (T       determine  each  other  uniquely  in  terras  of  the  Y 
•^    ^      '»  »"pti  pn  "^      '  pii 

which  we  may  assume  to  be  known.     The  initial  conditions  for  the    (f       are 

pji 

obtained  from  the  fact  that 


(9D) 


?  Cvr  •  [vf  ■  Cvr-    -^- 


As  before,  we  have 


Multiplying  by  any  R  ,  (50)  becomes 
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The    cr       are  predetermined  by  (23)  except  when  R     is  a  null  eigenvector  of 
H    *^  .     This  occurs  on  exactly  one  of  the  characteristic  surfaces,  say 
C^  .     Thus   (T       is  not  determined  by  (23)  and  its  initial  values  are  given 
by 

(53)  e''   (T*"    ♦  y''      -    <^''     -  r    (^eP(3^     ♦  y""  V  on   r. 

oqy'qv  ov^V^opv         pv  j  *  ' 

The  initial  values  of  the  cT"   will  be  smooth  functions  along  p  provided 


are  smooth,  i.e.,  provided  the  g    '  « 
guaranteed  by  our  assumption  that  g  is  smooth  on  both  sides  of  p ,  It 
follows  that  the  solutions  of  the  transport  equations  will  be  smooth  functions 
of  the  parameters  along  P,  hence  smooth  functions  of  position  on  the 
characteristic  surface,  (i.e.,  not  only  in  the  ray  direction). 

U.  Existence  and  uniqueness  theorem  for  the  discontinuous  initial  value  problem 

For  the  proof  of  the  theorem  which  we  are  about  to  present,  it  is 
necessary  that  we  be  able  to  solve  the  transport  equations  (plus  initial 
conditions)  along  the  characteristic  surfaces  C^  emanating  from  the  mani- 
fold p •  It  is  clear  from  the  foregoing  that  this  can  be  done  in  the  region 
0  <  t  <  t^,  where  t-  is  the  smallest  value  of  t  for  which  one  of  the  hyper- 
surfaces  C  '^  has  a  singularity.  Singularities  of  characteristic  surfaces 
bearing  non-zero  jumps  will  be  referred  to  as  caustics.  Whether  or  not 
caustics  appear  in  a  given  initial  value  problem  can  not  be  determined  solely 
by  examining  the  solutions  0  of  the  single  first  order  partial  differential 
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equation  (12)  whose  solution  sxa'faces  0*0  contain   p;  for  it  is  possible 
that  the  function  g(x)  may  be  such  that  although  the  characteristic  surface 
C^     emanating  from   P    has  a  singularity,  that  surface  does  not  play  a  role, 
the  corresponding  initial  values  for  the  transport  equations  being  all  zero* 
It  is  convenient  to  foiraulate  the  following 

caustic  condition:     The  problem  is  such  that  the  characteristic  surfaces 
across  which  u  and  its  derivatives  have  non-zero  Jumps  are  free  of  caustics. 
We  nay  now  state 

Theorem  2t     Given     the  initial  value  problem 

(5U)  Lu     -     0 

(55)  u(0,x)     -     g(x). 

Here  L  is  a  symmetric  hyperbolic  operator  with  smooth  coeff iciert  s 
which  satisfies  the  multiplicity  condition  and  g(x)  is  smooth  on 
either  side  of  a  smooth  (n-1)  dimensional  manifold    p  ^^  ^he  plane 
t  ■  0,  while  g  and  its  derivatives  have  jump  discontinuities  across 
P«     There  exists  one  and  only  one  weak  solution  u  in  the   class   A, 
for  t  in  the  interval  0<t<t|»       uis  smooth  for  all  x  and  t 
in  that  interval  except  on  the   characteristic  surfaces  emanating 
from   p •     Across  tbese  surfaces  u  and  its  derivatives  have  jump 
discontinuities  which  satisfy  the  transport  equations.     If  the 
caustic  condition  is  sAisfied  we  may  take  t,    ■  -foo. 

Proof;     Let  v(t,x)  be  an  auxiliary  function      which  is  smooth  everywhere  in 
■5" 


t,   is  the  smallest  value  of  t  for  which  a  caustic  appears. 

ind  I 
[3] 


The  device  of  introducing  this  function  and  much  of  the  foregoing  con- 


struction are  due  to  R.  Courant  and  P.  Lax' 
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the  region  0  <  t  <  t.  except  on  the  characteristic  sxirfaces  C  **  emanating 
from  p .  ▼  is  to  be  so  chosen  that  its  Junps  and  those  of  its  derivatives 
across  these  surfaces  are  equal  to  the  unique  solutions  of  the  transport 
equations  and  initial  conditions  for  u.  Then  the  function  T  ■  Lv  is  smooth, 
for  by  construction  the  function  Lu  and  all  of  its  derivatives  are  continu- 
ous across  the  C   •  Let  u  -  v  +  w  where  w  is  a  solution  of  the  inhomogeneous 
equation 

(56)  Lw  -  -  T  . 
The  Initial  conditions  on  w  are 

(57)  w(0,x)  -  g(x)  -  v(0,x). 

By  construction  of  v,  both  T  ■  Lv  and  w(0,x)  are  smooth.     Therefore  by 
Theorem  1,  the  solution  w  exists  and  is  smooth.     It  follows  that  u(0,x)   ■  g(x) 
and  Lu  ■  0  everywhere  except  on  the   surfaces  C  '^     where  u  satisfies  the 
transport  equations.     Now  as  in  (9)  and  (10), 

(58)  I    j  uM  i^dxdt  -   J   J  i\jLudxdt  +  J2     j       X1h°^P^  [u]PadS  . 

c 

But  Lu  ■  0,  and  H°'^   [u]^  •  o.     Hence  (U)  is  satisfied.     Thus  it  is  proved 
that  there  exists  a  weak  solution  u» 

For  the  uniqueness,  let  u  ■  u,   "  vu,  where  u,   and  Up  are  solutions  in 
^ ,    Then  u  is  in    ^   and  is  a  weak  solution  of  the  problem  with  g(x)   5  0» 
The  discontinuities  of  u  and  its  first  derivatives  can  occur  only  on  character- 
istic sxirfaces.     These  always  intersect  the  plane  t  ■  0.     But  the  initial 
values  for  the  transport  equations  on  these  surfaces  are  zero.     Hence  the 
discontinuities  are  zero.     It  follows  that  u  is  continuously  differentiable. 
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so  u  is  a  strong  solution.     Therefore,  from  the  imiqaeneas  assertion  of 
Theorem  1,  u  5  0.     This  proves  the  uniqueness  and  completes  the  proof  of 
the  theorem. 

5,     Initial-boundary  value  problems  for  symmetric  hyperbolic  equationa 
In  addition  to  initial  value  problems  for  syumetric  hyperbolic 
equations,  we  may  also  consider  the  so-called  mixed,  or  initial -boundary 
value  problems.     If  p(x)   "  0  is  the  equation  of  any  smooth  surface,  B,  in 
the  x-space,  the  transformation  z     »  a^Cx);  v  -  l,,,,,n;  with  o.     '  ^t  maps 
the  surface  B  into  the  plane  z     "0,  while  preserving  the  symmetric  hyper- 
bolic nature  of  L  and  the   smoothness  of  its  coefficients.     Thus  without 
loss  of  generality  we  may  take  our  boundary  to  be  x     "0.     On  this  boundary 
we  will  require  the  solution  u  to  satisfy  a  certain  type  of  homogeneous 
boundary  condition     which  we  will  now  explain. 

The  vector  u  is  an  element  of  a  k-dimensional  Fuclidean  vector 
space  V.     Consider  a  subspaoe  T  of  V  with  the  following  properties 

1.  The  matrix  a"  is  non-positive  over  T,  i.e.,  uA^u  <  0  for  every 
u  in  T. 

2.  T  is  maximal,  i.e.,  the  dL-nension  of  T  is  as  large  as  the  dimension 
of  any  subspace  having  property  1. 

The  boundary  condition  may  now  be  expressed  simply  as 

(65)  u(x,t)   in  T,  for        x^  -  0. 

If  T^  is  the  orthogonal  complement  of  T,  and  T     is  spanned  by  the  vectors 


In  Kaxwell's  eouations,  the  conditions  on  E  and  H  at  the  surface  of  a 
perfect  conductor  can  be  shown  to  be  conditions  of  this   type.     See  Part 
II,  Section  2, 
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b  ,*..,b  ,   the  condition  nay  also  be  stated  as 
(66)  ub  ■  Oj      8  •  l,,,,,ffl,   for  X    ^  0, 

8 

Here  ub     denotes  the  inner  product.     The  vectors  b  ,  and  the  space  T  need 

s 

not  be  fixed  but  may  be  functions  of  the  point  on  the  boundary.  We  will 
take  them  to  be  smooth  functions.  For  problems  involving  a  boundary  we 
will  make  use  of  the  following  generalization  of  Theorem  i:L°J»L'J  . 

Theorem  l-Bt  Let  L  be  a  symmetric  hyperbolic  operator  with 
smooth  coefficients,  let  (^(t,x)   be  a  smooth  function  defined 
in  the  domain  /^  j  t  >  0,  x  >  0;  and  g(x)  be  a  smooth  function 
defined  in  the  domain  x  >  0  which  satisfies  the  boundary  condi- 
tion: g  in  T,  for  x  -  0  and  t  -  0,  Then  the  initial-boundaiy 
value  problem 

Lu  -  0j   u(0,x)  -  g(x)}  u  in  T,  for  *  •  0 

has  a  unique  smooth  solution  u(t,x)  defined  in  ^  •  The  value 
of  the  solution  u  at  a  given  point  depends  on  the  values  of  the 
coefficients,  the  functions  0  and  g,  and  the  boundary  conditions 
only  in  a  finite  region. 

Weak  solutions  will  be  defined  as  before,  the  only  acl-ditional  require- 
ments being  the  boundary  condition  (6$)   and  the  requirement  that  the  region 
R  idilch  is  the  support  of  the  testing  function  SL.  shall  not  intersect  the 
half -space  x  <  0. 

The  entire  discussion  in  Sections  1,  2,  3  now  remains  valid  in  the 
i-egion  0  <  t  <  ninCt-jtp)  where  tp  is  the  smallest  value  of  t  for  which  one 
of  the  characteristic  surfaces  emanating  from  p  intersects  the  boundary. 
In  order  to  extend  our  existence  theorem  beyond  tj  we  mast  first  extend 
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our  a  priori  calculation  of  the  discontinuities  of  u.  This  is  done  in  the 
next  section, 

6,   Initial  conditions  for  the  transport  eqiaations  at  a  boundary 

Suppose  that  a  characteristic  surface  C  emanating  from  the  manifold 

P  in  the  initial  hyperplane  intersects  the  boundary  B  in  an  (n-1)  diaen- 

sional  manifold  ^.  /^  will  be  smooth  in  the  region  t  <  t, ,  for  there  C  is 

smooth.  We  now  consider  all  the  characteristic  surfaces  which  intersect  B 

in  ^  and  attempt  to  choose  initial  conditions  for  the  transport  equations 

on  these  surfaces  in  such  a  way  that  the  jumps  in  u  and  its  derivatives  will 

be  compatible  with  the  boundary  conditions  (65). 

Let  5  ,,..,5  ,,0  denote  components  of  the  unit  normal  to  /^  in  the 
o    '  n  ~x 

plane  x  •  0.  Let  0°^^'  ■  0  be  the  equation  of  a  characteristic  surface 
n 

C^^^  whose  intersection  with  B  is  A.  Then  <;S°^^'   ■  pfi^,  on  A  (v  -  0,...,n-l) 

where  p  is  a  proportionality  factor.  Here     2S^^^  is  a  solution  of 

n 

(67)  det  H°^P^  -  det/^n    ^x^^^  ^'']    "  ^et  (V^^^^  a"  *  pf)  -     0}  on  A 

\v»o       V  /  \    r*  / 

where 

n-1 


(68)  F    -    n    S^  A 

v»o 


V 


Equation  (6?)  will,   in  general,  have  complex  roots  0      ^^  ,     For  each  distinct 

n 
real  root  there  will  be  a  corresponding  characteristic  surface  through  A* 

We  assume  that  the  manifold    p  is  bounded  away  from  the  boundary,  B.     Then  it 

is  easily  seen  that  there  will  be  a  point  P  on  A  °^  minimum  t,   and  that  at 

p,  5,   .....  6     ,   -  0.     Hence,   at  P,  6     ■  1  and  F  ■  I.     It  follows  from  the 
*     1  n-1  o 
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continuity  of  Fj  that  in  some  neighborhood  of  P,  F  is  non -singular.  We  now 
Inpose  a  new  condition  on  the  problem}  the 

shadow  condition:  F  shall  be  non-singular  throughout  ^«  This  condition 
as  suggested  by  its  name  has  an  interesting  geonetric  interpretation  which 
will  be  explained  below.  It  follows  fron  the  shadow  condition  that  F  is 

positive  definite  on  ^.  Iherefore  we  may  set  F  ■  JJ  where  J  is  non- 

n  ' 
singular  and  has  an  inverse  S»  Let  G  -  SA  S  •  Then  G  is  syimnetric  and 

we  may  write 

(69)  (G-ttp)  Q^P^  -  Oj    p  -  l,...,k. 

Here  the  q''^'  are  a  complete  set  of  prthogonal  eigenvectors  of  G,  which  we 
may  normalize  by  the  condition 

(70)  q^P^Ss'q^P^  -  Ij         p  -  l,..,,k. 

The  eigenvalues  a  are  real.  Set  R^^'  -  S  Q^^'j  p  •  l,..»>k.  It  follows 
from  (70)  that  the  R^^'  are  unit  vectors.  They  are  clearly  linearly  inde- 
pendent. From  (69),  J(SA^s'-a  )j'r^^^  -  0,  or 

(71)  (A'^-«pF)R^P^  -  0}   p  .  l,...,k. 

For  each  non-zero  a  ,  set  9^  **  ■  "P/a_>  then 

P       n        ^ 

(72)  H°^P^R^P^  .  C(^°^P^  a"  *  p?)   R^P^  -  0}   if  a^  jt  0. 


Equation  (72)  expresses  the  fact  that  for  each  distinct  non-aero  a  ,  there  is 
a  distinct  characteristic  surface,  C^^  ,  passing  through  A»  and  R^^  is  a 
null  eigenvector  of  the  associated  matrix  H°  ^  .  Since  we  may  have  multiple 
roots  a  ,  we  may  have  several  distinct  R  ^  associated  with  the  same  matrix 
H°^P  ,  As  we  have  seen  earlier,  the  jump  in  u  across  C^^  ,  [u]^,  can  be 
expressed  in  terms  of  the  associated  null  eigenvectors.  Therefor* 
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(73)  [u]  ^  -  E  w  "*  ■  i:  <  R^^^  • 

q  P-1     ^ 

Since  we  have  used  only  the  non-zero  o  ,  vre  may  assume  that     <f  °  •  0  in  (73) 
if  o     ■  0.  (if  A     is  singular  there  will  be  zero  eigenvalues  a   .     For  these 
o    we  have  solutions  of  (72)  with  p  ■  0,     In  this  case  the  boundary  is  itself 
a  characteristic  surface.)  We  apply  the  boundary  condition  in  the  form  (66): 

A        .Vr.A         k. 
P" 


(7U)  [b^]^.  b^[u]^.  f    bVPV°  .0}  V  .  l,...,m 


Of  course,  some  of  the  (5^  are  already  determined  (roughly  speaking,  those 
associated  with  a  characteristic  surface  emanating  from  the  initial  plane 
or  from  the  boundary  at  an  earlier  time).  Our  object  will  be  to  show  that 
(7U)  can  be  solved  uniquely  for  the  undetermined  (f  in  terms  of  the  deter- 


mined ones.  It  is  easy  to  see  that  (T 

P 


•"will 


will  not}  ^^""^   ®  preas signed 
value  at  a  point  on  Z^  if  the  ray  at  that  point  (oriented  in  the  direction 
of  increasing  t)  proceeds  ■°y*'.°^[  the  region  x  >  0,  i.e.,  if  x  -  R^P'a"r^P' 


is  ■ 


^  Q  \    •     (V!e  shall  see  that  the  possibility  x     ■  0  is  precluded  by  the 

shadow  condition.) 

We  digress  to  explain  the   significance  of  the  shadow  condition: 

Since  the  R  *^     are  linearly  independent  they  maQr  be  made  the   columns  of  a 

non-singular  matrix  R,     Then  the  matrix  R  FR  »  M  -  (m     )  has  entries 

»pq  -  R^pW'I^     But  R^pW^^^  -  R^P^Tt'r^^^   -  Q^pV^   -  0  if  p  ^  q. 

Hence  M  is  a  diagonal  matrix  and  det  M  ■  "ff    R  P  FR^P\     The  shadow 

p-1 
condition  now  guarantees  that 

(75)  R^pW^P^     ^0}  p  -  l,...,k  . 
From  (71), 

(76)  r(^^A°R^P^    •-     aR^^W^P^ 

P 
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For  all  characteristic  surfaces  C^^     except  the  boundary  itself  a     /  0, 
hence  x^  -  R^^^A^R^^^  /  0.     But  i^  ■  0  if  and  only  if  the  ray  is  tangent 
to  the  boundary.     Thus  the  shadow  condition  suffices  to  ensure  that  no 
ray  will  be  tangent  to  the  boundary.     In  the  language  of  physical  applic- 
ations, we  say  that  there  are  no  shadows  in  the  problem. 

Returning  to  our  examination  of  the  linear  algebraic  system,   (7U), 
let  RfPU-^'P)  be   {>  0  f-  P  :  J:i;:;:,,)    .    we  r,.rit,  (7U)  placing 
known  quantities  on  the  right  side 

(77)  t    bVPV°     .  -  JZ      bVPV°  ;       V  .  l,...,ra  . 

P«l  P       paS+1  " 

Our  object  of  proving  that  (77)  has  a  unique  solution  will  be  attained  if 
we  can  show  that  s  ■  m  and  the  coefficient  matrix  (b  R  ^  )  is  non-singular. 
To  this  end,  we  shall  employ  some  simple  lemmas  proved  in  Appendix  2.  From 

(76),  R^pW^^^  .  0  if  p  /  q.  Let  V  .  f  c  R^^\     Then  va"v  -  f^  c^^P^a"r^P^>  0 

(1)      (s)      P=l  ^  P^l  ^ 

unless  V  ■  0,  Hence  R   ,...,R    span  a  subspace  T,  over  which  A  is  positive. 

Similarly  R     ,...,R    span  a  subspace  T^  over  which  A  is  non-positive. 

It  follows  from  lemmas  U  and  2  that  dim  T,  <  m  and  dim  T«  <  k-m  where  m  is 

the  number  of  positive  eigenvalues  of  A  .  From  lemma  3,  m  is  also  the 

dimension  of  T   (which  is  as  it  is  used  in  (66),  (7U),  and  (77)).  But 

dim  T,  +  dim  Tp  ■  k.  Hence  dim  T,  ■  s  •  m.  Thus  the  coefficient  matrix 

(b^R^P  )  of  (77)  is  square.  Suppose  it  were  singular.  Then  there  would 

m  I  (    ^ 

exist  non-trivial  numbers  d     such  that  T"    d^  b  R^^^  «  0;  v«l,..,,m.     Set 

p         *— 1  p 

u  -  7"  d  R^P\  Then  ub^  -  Oj  v-l,...,mi  and  ua"u  -  T"  d^R^P^A^R^P^  >  0. 

p^  P  P^  ^ 

But  this  contradicts  the  boundary  conditions  (see  Section  5)«  Hence  the 

coefficient  matrix  is  non-singular  and  the  solution  of  (77)  exists  and  is 

unique.  That  solution^  <5'^,..»>  (fZ.   »  provides  the  initial  values  for  the 

transport  equation  of  order  zero. 
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To  obtain  the  initial  conditions  for  the  higher  transport  equations 
we  proceed  as  followss     Let  s  be  a  unit  vector  in  the  ray  direction,   and 
5-  denote  differentiation  in  that  direction.     Now  s  and  y     being  orthogonal 
directions,  they  determine  a  two-dimensional  plane  which  must  intersect  the 
(n-dinensional)  boundary  in  a  direction,  N  which  lies  in  the  boundary,  but 
which  cannot  lie  in  zu»     This  is  true  because  y     is  orthogonal  to  Z^  and  the 
shadow  condition  guarantees  that  s  is  not  in  ^,  (see  Fig,  2).     We  now 
differentiate  the  boundary  condition  (66)  in  the  direction  N,  r  tines  to 
obtain 

(7B)  u      b^  +  ru  /     -  N  b„  +   •••     «     Oj     on  x     -0;     v  «  l,.,,,in} 

jjT  N^         -^ 


or 

(79)  ju  J       b^  -  -rju  (j._i)J       bJJ  -  •••;  on  x^  -  0;     v  -l,...,m. 


n 


As  in  Section  2,  we  express  u       in  terms  of  derivatives  in  the  y     and  s 

N 
directions: 

r 

(80)     u    .  7"   e^  u  ,   ^ 

„r    ^—    V   (r-v)  V 
N      v-o     y^    's 
o 


or 


'«^'    [vJ'-'UvT'fe'v  [»<-,] 


y^-,    v-1    ^  y^   '-.s 


It  is  important  to  note  that  the  shadow  condition  guarantees  that  the 
directions s  and  N  do  not  coincide,  hence  e^  /  0.  Suppose  now  that 


•'o 


o 

p  (p) 

has  been  determined  for  each  characteristic  surface  C    which 


intersects  the  boundary  in  the  manifold  Zi^,  for  p  -  0,1,..,, r-l.  Then 
the  jumps  in  al3  derivatives  of  order  less  than  r  are  known.  To  avoid 
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Figure     2 


n 


Relative  orientstions  of  the  directions  y  ,  N,  s  at  the  boundary. 
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unnecessary  writing,  let  w  denote  known  quantities.  Then 


'-   [v]^-frvT-F'°[vr 


From  (79) 

iP 


(83) 


n    eM   u  ^         b^     -    ^2;  V  -  l,...,m. 


On  each  C       ,       u  can  be  expressed  in  terms  of  the  eigenvectors  of 


o(   )  ^° 

H    "  ,  the  coefficients  of  the  eigenvectors  which  are  not  null  being 

known  from   (23).     The  totality  (for  all  C^^O  of  the  null  eigenvectors 


t^P^   p  »  l,...,k. 


are  simply  the  R  ^  j   p  *  l,,,,,k.     Furthennore  the   coefficients  of 
R  ,,.,,R         are  knownj  hence 


(8u)  Yl    ^^  ^^"^^  ^l  K  '  ^y     V  =  1,...,™  . 

p-1  p 

Comparison  with  (77)  now  shows  that  (8U)  can  be  solved  uniquely  for  the 
quantities  e^  q   ;  p  -  l,...,ni;  hence,  since  e  /  0,  for  the  <^   ,  which 
provide  the  initial  values  for  the  transport  equation  of  order  r. 


7,  Existence  and  uniqueness  theorem  for  the  discontinuous  initial  boundary 
value  problem 

In  order  to  generalize  the  theorem  of  Section  3  to  problems  with 
boundaries  we  need  only  impose  upon  the  auxiliary  function  v  the  additional 
requirement  that  it  satisfy  the  boundary  condition.  This  can  be  done  because 
the  discontinuities  across  the  characteristic  surfaces  which  intersect  the 
boundary  are  such  as  to  be  compatible  with  the  boundary  conditions.  Now, 
by  Theorem  1-B  we  can  require  that  the  function  w  also  satisfy  the  boundary 
condition.  Then  the  function  u  »  v  ♦  w  will  satisfy  the  boundary  condition. 
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Ihe  uniqueness  argument  requires  an  obvious  modification*  Thus  we  have 
Theorem  2-Bi  Given  the  initial  boundary  value  problem 

(85)  Lu  -  0 

(86)  u(0,x)  -  g(x) 

(87)  u(x,t)  in  T,    for  x   -  0. 

n 

Here  L  is  a  syimnetric  hyperbolic  operator  with  smooth  coefficients 

defined  in  the  domain  /0:t>0;x  >0;  g(x)  is  smooth  on  either 

—    n  — 

side  of  a  smooth  (n-1) -dimensional  manifold  p  ^  ^^®  plane  t  ■  0, 
while  g  and  its  derivatives  have  Jump  discontinuities  across  P 
and  g  satisfies  the  boundary  condition:  g  in  T  for  x  -0  and 
t  ■  0.  Let  the  multiplicity  and  shadow  conditions  be  satisfied. 
Then  there  exists  one  and  only  one  solution  u  in  the  class  A 
in  the  region  0<t<t.,x  >0.  uis  smooth  everywhere  in 
that  region  except  on  the  characteristic  surfaces  emanating  from 
\     and  from  the  intersections  i\   (which  may  be  multiple)  of  such 
surfaces  with  the  boundary,  x  »  0,  Across  these  surfaces,  u  and 
its  derivatives  suffer  jump  discontinuities  which  satisfy  the  trans- 
port equations.  If  the  caustic  condition  is  satisfied  we  may  take 
t^  ■  ♦00, 
note:  By  virtue  of  the  remarks  made  at  the  beginning  of  Section  5  the  above 
theorem  is  immediately  valid  for  any  domain  aj    which  can  be  mapped  by  a  smooth 


« 
It  is  clear  that  our  construction,  and  hence  this  theorem  are  vaHd  for 


multiple  'reflections'  from  the  boundary. 
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* 

transfonnation  of  the  space  coordinates  onto  the  domain  t>0;x    >0, 

—     '     n  — 

The  theorem  is  also  valid  for  more  general  domains  such  as  the  interior 
of  a  cylinder  with  generators  parallel  to  the  t-aixis  or  domains  botinded 
by  several  smooth  surfaces.     In  these  cases  we  must  employ  the  analogues 
of  Theorem  1-B  valid  for  such  domains'-  J'L'J^ 


7 

It  should  be  noted  that  the  mapping  induces  a  change  in  the  statement  of 

the  boundary  condition  and  the  shadow  condition*  Thus,  e.g.,  if  p(x)  »  0 

defines  the  boundary^the  boundary  condition  becomes:  u(x,t)  in  T  for 

n 
p(x)  «  0,  where  T  is  the  maximal  subspace  such  that  the  matrix  7~  p  A 

J.S  non-positive  over  T. 
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Appendix  1 

Lemma;     Let  P,  A  be  real  syinraetric  matrices  of  dimension  kj 
let  M(^)   =  P  +  j^A  have  eigenvali»s  whose  multiplicities  are 
the  same  for  every  real  value  of  ^,     Let  X(^)  be  an  eigen- 
value of  M  with  associated  orthonormal  eigenvectors  R  (^), 
,..,r''(^).     (Thus  q  is  independent  of  >^,)     Then  X(^)  is  an 
analytic  function  of  }^  in  a  neighborhood  of  every  real  point 
Y.f  and  the  R  (^)  may  be  so   chosen  as  to  have  the  same  pro- 
perty. 

Proof;     \  is  a  solution  of  the  equation 

(1)  det(XI-M)  .  G(\,^)   -  X^  +  e(k-l)^^^^^'^*   •*•  *  ^0^^^  "  °* 
Here  the  g     are  polynomials  in  S^.     Let  us  factor  G  into  irreducible  poly- 
nomials H  of  the  form 

(2)  H(X,>:)  -  hp(^)xP  +  %.i)(^)^^*^*  •••  ♦  h^(^)« 
Then  each  eigenvalue  X  is  a  solution  of  an  equation 

(3)  H(X,^)     -     0. 

Thus  X  is  by  definition  an  algebraic  function  of  ^,  From  the  theory  of 
algebraic  functions  ,  we  know  that  the  equation  H(X,^  )  ■  0  has  p  distinct 
roots  except  at  certain  exceptional  points  5^  where  either 

1,  h  vanishes 

P 

or 

2,  the  discriminant  of  the  equation, which  is  a  polynomial 
in  J^,  vanishes. 

The  discriminant  of  an  irreducible  polynomial  cannot  vanish  identically. 

Hence  the  singular  points  of  type  2  are  finite  in  number.  Furthermore  in 

See,  e.g.,  [8^  for  the  statements  about  alpebraic  functions  which  we 

employ  here. 
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our  case  there  can  be  no  real  points  of  type  2,  for  at  such  points  the 
multiplicity  of  some  eigenvalue  must  be  greater  than  at  neighboring 
points.  Also  we  observe  that  h  (i^)  must  be  a  constant,  so  there  can 
be  no  points  of  type  1,  Since  algebraic  functions  are  analytic  except 
at  exceptional  points,  the  assertion  for  X(^)  follows. 

The  eigenvectors  R  are  solutions  of  a  linear  algebraic  system 

(U)  TR  -  0. 

Here  T  ■  (^^J  "  MK)1   -  P  -  ^A,  Thus  the  t.  ,  are  analytic  functions  of 

>^  for  real  ^.  The  nullity  of  T  is  q,  and  if  r  denotes  the  rank  of  T,  we 

have 

(5)  r  +  q  »  k. 

To  solve  (U)  for  the   components  of  R  we  may  first  reduce  it  to  an 
equivalent  system  of  r  equations  and  k  unknowns.     We  then  assign  arbitrary 
values  to  k-r  -  q  of  the  unknowns  x, ,..»,x     such  that  the  remaining  r  x  r 
coefficient  matrix  is  non-singular.     We  may  then  solve  for  the  remaining 
r  unknowns  by  Cramer's  rule.     The  x     ^,...,x.    will  therefore  be  rational 
functions  of  the  x^,...,x     and  the  t .  .  with  non-vanishing  denominator.     By 
choosing  x.  ,...,x     as  analytic  functions  of  J^  we  are  guaranteed  that  R  is 
an  analytic  function  of  ^,     We  obtain  q  linearly  independent  solutions 
B.iK)t  snd  the  analyticity  is  preserved  in  the  orthonormalization  procedure 
to  obtain  the  R^(^),.  ..,R^(^). 
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Definitions:     Given  a  k-dimensional  vector  space  V,  a  real  syinmetric  operator 
A,   and  a  sub  space   T.     A  is  non -positive   over   T  if  (vi^u)  <  0  for  all  u  in  T. 
A  is  positive  over  T  if  (u,Au)  >  0  for  all  u  /  0     in  T,     Let  a  ,  r     denote 

V  V 

(real)  eigenvalues  and  orthonormal  eigenvectors  of  A;  i.e.,  Ar     »  a  r  , 
V  ■  l,.«.,k.     Let  a     >  0  for  v  ■  l,,..,m;  a     <  0  for  v  =■  m+l, . .  .,k, 

Leitima  1;     Let  T,  R  be  disjoint  subspaces  of  V  (i.e.,  if  u  is  in 
T  and  in  R  then  u  •  0),  then  dim  T  +  dim  R  <  k. 

Proof;     The  proof  of  this  lemma  is  trivial. 

Lemma  2;     Let  A  be  non-positive  over  T.     Then  dim  T  <  k-m. 

Proof!     Let  R  be  the  subspace  spanned  by  r^,,,.,r   •     Then  T  and 

R  are  disjoint,   for  if  u  is  in  R  and  in  T,  we  have 

m 
"     -    II    «v^  ' 

A'^     -    II    '^v^v'"^' 

2 
(u,Au)     -     l_    c^a^     <     0. 

Hence  c     ■  0;   v  «  l,,..,m.     Therefore  u  ■  0, 

Lemma  3;     Let  T  be  a  subspace  of  V  of  maximal  dimension  such 
that  A  is  non-positive  over  T.     Then  dim  T  =  k-m. 

Proof;     Let  T,  be  the  subspace   spanned  by  %+]>•• '^^k*     T^^^^  A 

is  non-positive  over  T,   and  dim  T,   -  k-m. 

Lemma  U;     Let  A  be  positive  over  P.     Then  dim  P  <  m. 

Proof:     Let  T,  be  defined  as  above.     Then  T^  and  P  are  disjoint, 

for  if  u  ^  0  is  in  T,    and  in  P,     we  have 

"     ■     II  ,      ^v^   i 

v=m+l 

(u,Au)     -     5Z     c^a^     >     0, 


which  is  a  contradiction. 
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Part  II 

Asymptotic  expansion  of  steady-state  solutions  of  symmetric  hyperbolic  linear 
differential  equations 

Introduction 

fol 
In  a  series  of  lectures'-  J  given  at  New  York  University  during  19U7- 

U8,  R.  K.  Luneburg  introduced  a  method  for  obtaining  asymptotic  expansions 
of  solutions  of  Maxwell's  equations.  The  theory  was  refined  and  extended  to 
linear  second  order  hyperbolic  partial  differential  equations  by  M,  Kline L  J. 
Similar  asymptotic  expansions  have  been  studied  by  several  authors 
and  have  been  applied  to  the  solution  of  diffraction  problems L  -1.  These 
methods  have  the  virtue  of  providing  the  asymptotic  expansion  of  the  solution 
even  where  it  is  difficult  or  impossible  to  obtain  the  solution  itself.  In 
cases  where  the  exact  solution  is  known,  the  expansion  is  frequently  more  use- 
ful and  can  be  obtained  more  easily  by  these  methods  than  by  expanding  the 
known  solution. 

The  present  paper  is  concerned  with  extending  the  Luneburg -Kline  theory 
in  two  directions:  First,  we  will  consider  more  general  partial  differential 
equations.  Second,  we  will  attempt  to  provide  proofs  of  certain  assertions  for 
which  Luneburg  and  Kline  use  plausibility  arguments.  Our  goal  is  a  rigorous 
proof  of  the  asymptotic  nature  of  the  formal  expansions  which  will  be  derived. 
In  this  proof  the  results  of  Part  I  play  a  key  role  but  nevertheless  we  have 
attempted  to  make  the  present  discussion  self-contained. 

The  plan  of  presentation  is  as  follows:  In  the  remainder  of  this  in- 
troduction we  discuss  notation  and  the  formulation  of  the  problem.  Then  in 
Section  1  a  formal  expansion  is  derived.  In  Section  2  this  expansion  is 
applied  to  Maxwell's  equations.  Thus  Sections  1  and  2  contain  our  formal 
results.  The  remaining  sections  are  concerned  with  the  proof  of  the  asymptotic 
nature  of  these  formal  results. 

We  now  turn  to  the  definitions  and  formulation  of  the  problem.  A  fuller 


See  the  references  listed  in  Qfj  ■ 
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discussion  of  some  of  the  following  notions  is  contained  in  Part  I.  We  denote 

by  u, a  k  component  vector  function  of  the  (n+1)  variables  x  ,  x.,...,x  ,  and 

by  A  and  B, k  x  k  matrices  with  real  entries.  We  will  sometimes  denote  x  by 

t  and  x.  ,...,x  by  x.  The  matrices  A  and  B  are  taken  to  be  smooth  fxinctions 
1'   '  n  ' 

of  X  in  the  entire  x-space.  We  define  a  first  order  linear  differential 
operator 

(1)  Lu  -  |~  a\       +  Bu        (u   •^). 

v-o    *v  V   "^-^v 

V 

The  operator  L  is  taken  to  be  symmetric  typerbolic,  i.e.,  the  matrices  A  are 
assumed  to  be  symmetric,  and  the  matrix  A  to  be  positive-definite.  Without 
loss  of  generality  we  may  take  A  •  I.  Now  we  consider  the  initial  value 
problem  for  the  vector  w(t,x) 

(2)  Lw  -  g(x)e"^"'*  }  w(0,x)  -  h(x). 

Here  g  and  h  are  vector  functions  and  u  is  a  (scalar)  parameter  which  will 
be  referred  to  as  the  (circular)  frequency.  In  physical  applications,  the 
function  ge~    represents  a  source  with  harmonic  time  dependence. 

Under  appropriate  conditions  (which  will  be  made  explicit  later)  the 
solution  w  of  (2)  may  be  decomposed  uniquely  into  two  terms: 

(3)  w(t,x)  -  z(t,x)  +  e"^^  v(x) 

where  z(t,x)  ->  0  as  t  ->oo.  z  is  called  the  transient,  and  v(x),  will  be 
referred  to  as  the  steady-state  function.  The  problem  we  consider  is  that 
of  obtaining  the  asymptotic  expansion  of  v  with  respect  to  w,  as  oo  ->oo. 


The  term   'smooth'  will  be  used  to  describe  functions  which  have   continuous 
derivatives  of  sufficiently  high  order  for  the  purpose  at  hand. 
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It  will  be  shown  that  v(x)  is  independent  of  the  initial  I'anction  h  and  is 
a  solution  of  the  reduced  equation j 

n 

(U)  Nv  »■  ^  A  V   +  (B-ico)v  ■  g, 

v-1    ^v 

(Fq.  (U)  may  be  obtained  formally  by  substituting  w  ■  e    v  into  eq.  (2).) 

Although  the  reduced  equation  does  not  suffice  to  determine  v  uniquely,  it 

can  be  used  in  formally  deducing  the  expansion  of  v.  In  this  deduction  it 

is  necessary  to  impose  on  the  expansion  an  outgoing  condition  which  will 

be  explained  later.  In  the  later  sections  it  will  be  shown  that  the  formal 

expansion  is  indeed  an  asymptotic  expansion  (for  co  ->oo)  of  v(x)« 
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1,  Ttie  formal  expansion 

VJe  seek  a  fonnal  expansion  v  of  the  steady  state  function  v  defined  in 
the  Introduction,  This  function  satisfies  the  equation 

(5)  Nv  -  g(x) 

Our  procedure  will  be  formal}  the  results  will  be  justified  later.  We  begin 
by  considering  an  expansion  of  the  form 

(6)  y  ■  E    ^  . 

r»o  ( ico; 

(r) 
Here  the  y       (x)   are  vector  functions  to  be  detemdned.     If  we  insert  (6)   in 

(5)  and  compere  powers  of  co,  we  obtain 

(7)    ■  y^°^-Oi         y^^^--g}         y^"*^^-/y^"^         r  -  1,2,...     . 

The  operator  oc   is  given  by 

(8)  /  .  r  A^  ^ 


♦  B. 


^\ 


(r) 
From  (7)  we  obtain  the  functions  y^  ^   explicitly  in  the  form 

(9)        /°'.0}     /""^  '  -  ^^'-^\i  r.1,2,...  . 

Here  ^°  -  I.  Thus 

00    y>r 

(10)  y   -r    7^. 

r»o  (ico; 
If  g(x)  is  infinitely  differentiable,  the  desired  formal  solution  is  simply 

(11)  V  -  y. 

If  the  source  region  is  bounded  (i.e.,  if  g  5  0  outside  a  bounded  region)  the 
formal  solution  vanishes  identically  outside  that  region. 
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Let  us  now  consider  the  more  interesting  case  in  which  g  or  one  of  its 
derivatives  is  discontinuous.  More  precisely,  let  p  be  a  smooth  (n-1)  dimen- 
sional manifold  in  x-space,  and  let  g(x)  be  smooth  on  both  sides  of  P ,  while 
g  and  its  derivative?  have  jump  discontinuities  across  f  ,  Later,  in  consider- 
ing problems  with  boundaries,  we  shall  require  that  g(x)  vanish  in  a  neighbor- 
hood of  the  boundary,  A  function  g(x)  with  all  these  properties  will  be  said 
to  be  of  type  g.  We  require  that  our  formal  solution  v  shall  be  continuous, 
(It  will  be  seen  later  that  the  exact  solution  ▼  is  also  continuous.)  To 
achieve  this  we  may  begin  with  (10)  but  we  must  repair  the  discontinuities 
in  it  without  disturbing  the  equation  (5).  We  nay  hope  to  accomplish  this 
by  adding  to  y  solutions  of  the  homogeneous  equation 

(12)  Nv  -  0. 

For  this  purpose  we  shall  study  functions  of  the  form 

(13)  z^  =  ^  -  e"'^   >''■'  JZ      5 _iH  ^ 

r-o    (ico) 

Here  a  is  an  index,  y'°'  Is  a  scalar  function,  z^°  and  2^°"'^^     are  vector 
functions  . 

If  we  apply  the  operator  N  to  (13)  and  equate  the  coefficient  of  each 
power  of  (io))  to  zero  we  obtain 

(lU)  G^^^  2^°^^°^  -  0 

(15)  G^°^  ,(a)(r.l)  .  /,(a)(r)  .  q.        ^  ^  0,1,2,... 

where 

The  formal  expansion  y  defined  by  (6)  is  a  special  case  of  the  form  (13).  In 
fact,  if  one  attempts  to  satisfy  the  inhomogeneous  equation  (5)  by  inserting  (13) 
one  is  forced  to  set  Y    S  0.  Thus,  our  formal  procedure  requires  the  single 
''ansatz"(13). 
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(16)  G^"^  -  I  -  r  7^^^  A^  . 

If  the  matrix  0^°'^  were  non-singular,  these  equations  would  have  only  the 
trivial  solution  z^°^    '   5  Oj  r  ■  0,1,,.,     ,     Hence  we  assume  that 

(17)  det  G^°^  -  0, 

Before  taking  a  closer  look  at  the  equations  (lU),  (15),  it  is 
necessary  to  impose  an  Important  condition  on  the  class  of  operators  L 
(or  N)  under  consideration,  which  we  shall  call  the 

/  \  r—  V 

multiplicity  condition!  If  P(x,z)  ■  ^  z  A  ,  where  the  z  are  real  numbers 
and  not  all  zero,  then  the  multiplicities  of  the  eigenvalues  of  the  matrix  P 
are  independent  of  x  and  z  «  z^,,..,z   , 

For  convenience,  we  shall  temporarily  omit  the  index  (a)  in  our  con- 
sideration of  a  particular  z       ,     From  (17)  we  see   that  the  matrix  G  is  singular, 
so  it  has  at  least  one  zero  eigenvalue  at  every  point.     At  some  point,  let 
X.,...,\     denote  the  zero  eigenvalues  of  0  (i.e.,  \   '  \  "  •*•  •  X     ■  0), 

Then  the  multiplicity  condition  assures  us  that  the  eigenvalues  X,(x),,,.,\  (x) 

*  la 

of  G(x)  are   zero  everywhere,  and  all  the  otters  are  non-zero   «     Let  R  ,«,,,R^ 


be  an  orthonormal  set  of  null  eigenvectors  of  G,  i,e.,  those  eigenvectors  out 

1  k 

of  the  complete  set  R  ,,.o,R     of  orthonormal  eigenvectors  which  are  associatec 

with  the  eigenvalue  X.   a  •••  ■  X     s  0,     From  (lU)  we  see  that  z^       is  in  the 
null-space  of  G,   so  it  has  an  expansion  in  terms  of  the  null  eigenvectors; 

(18)  z^°'     -    f    ^"r'" 

*— ,      m 

m»l 

(tVie  C^  «=  <r°(x)   are   scalars).     We  obtain  a  differential  equation  for  z^*^^ 
from  (1^')   for  r  =  0  by  multiplying  (15)  by  a  null   eigenvector: 

(19)  t     R^aV"^   -R^Bz^°^     -     Gj  i.l,...,q. 


v«l  '^v 


For  e  clarification  of  this  assert,ion,  see  Part  I,  Section  2,  the  second 
footnote. 
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Inserting  (l8)  in  (19)  we  obtain 
(20)   f  f  R^AV(<r°)^  .  f  f  R-^aV  ((r°)  .  f  R^  BR»<;  -  Oj  i-  l,...,q. 

Equation  (20)  is  a  first  order  system  of  linear  differential  equations  for  the 
scalar  coefficients  C    in  the  expansion  of  z^   ,  Before  discussing  it  further, 

we  will  derive  the  corresponding  equations  for  z   ,  r  >  0. 

(r) 
To  do  this  we  expand  z^   in  terms  of  all  the  eigenvectors  of  G 

(21)  z^^^  -  r   ff^  r". 

in-l 
From  (15)   and  (21)  we  have 


.(r)  .  ^    ,v     (r-1)  ^  ,  -(r-l)       ^ 


(22)  Gz^^^   -  T    a"  z^^  ■"'  *  Bz^^  -^^  .  Z.       o-m  V  }         '  "  ^*2,... 

v^  V  m-q+1 

Here  the  eigenvalue  X  is  simply 

(23)  X   -  r"  GR"". 

m 

Multiplying  (22)  be  an  eigenvector  which  is  not  null,  we  obtain 


(2U)        <r!:-i- Jr  RVz^r'^.RV^-^^ 

v«l 


«   \  1  t=l      ^v 


s  ■  q+l,.,.,k 
>  *         r  -  1,2,... 


Thus  we  see  that  the  coefficients  of  the  non-null  eigenvectors  in  the  expansion 
(21)  of  z"^'   are  given  explicitly  by  (2U)  in  terras  of  z^''"'^'o 

To  obtain  the  coefficients  of  the  null  eigenvectors,  we  multiply  (l5) 
by  a  null  eigenvector 

Inserting  (21)  in  (25)  and  rearranging  terms  we  obtain 
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(26) 


where 


X  ■  1, ...,q 


(27)     f,   -  r  r^aV  +  r^br"" 


and 

(28)  h    "t     t     R^AV(<r^)   -X:   ^s  <  • 

By  virtue  of  (2h),  g.  is  piven  explicitly  in  terms  pf  z     .  If  we  establish 
the  convention  z^    «  0  (20)  appears  as  a  special  case  of  (26).  (This  has  been 
anticipated  by  including  r  •  0  in  (26).) 

It  is  possible  to  put  (26)  in  a  much  simpler  form.  For  this  purpose  we 
employ  the  following  lemma  which  is  a  consequence  of  the  multiplicity  condition  . 

Lemma;  Let  R  ,,,.,R^  be  eigenvectors  associated  with  the  same 

n     ^ 
eigenvalue  \  of  the  matrix  P  ■  7°"  z  A  (z  .  real).  Then  for 

O  ^"'t       V      V 

v»l 
each  J, 

(29)  R'^  aV  -  c.  5  , 

o     o     J  Xm 

(i.e.,  R^A'^r"  -  0  if  /  ;^  m,  and  is  independent  of  /  if  /-  m.) 
It  follows  immediately  from  the  lemma  that  (26)  may  be  written 

(30)         ^ "W'-D^/ 1,  ^i.  <  •  4      tilling..  . 

We  now  introduce  a  family  of  curves  in  the  x-space  defined  parametrically  by 
x^     »     x^is),  where 


* 
The  proof  appears  in  Part  I,  Section  2, 
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(31)  ^v  '  "3r    ■  R  aVj  V  -  l,...,n. 

These  curves  will  be  referred  to  simply  as  rays«     Thenequation  (30)  may  be 
written  as  an  ordinary  differential  equation  along  a  ray: 

t  '  l,...,q 


(32)  ^,      ■  /_      -,7-.  -  «  -  6.  r»  0,1,2,... 


For  each  fixed  r,  (32)  is  a  system  of  q  linear  ordinary  differential  equationSo 
The  recursive  system  of  differential  equations  (32),  together  with  the  algebraic 
equations  (2U)  will  be  called  the  transport  equations. 

The  function  T  appearing  in  the  expansion  (13)  of  z    is  called  the  phase 
function.  We  have  seen  that  it  is  a  solution  of  the  single  first  order  non- 
linear partial  differential  equation  (17) »  a  special  case  of  which  is  called 
the  eiconal  equation.  It  is  interesting  to  note,  and  not  difficult  to  verify 
that  our  'rays'  are  the  'characteristic  curves'  used  in  constructing  the 
solution  of  the  first  order  equation  for  T.  We  also  note  that 

n  n 


(33) 


V^l    V       V«l      V 


Thus  T  may  be  identified  (up  to  an  additive  constant  on  each  ray)  with  the 
curve  parameter  s. 

A  formal  series  z    of  the  form  (13) »  where  T    is  a  solution  of  the 
eiconal  equation,  and  the  z      (or  rather  their  scalar  components)  are  solu- 
tions of  the  transport  equations  will  be  called  a  wave.  A  wave  is  not  completely 
determined  until  we  specify  initial  conditions  for  the  eiconal  e quation  and  for 
the  transport  equations.  Thus  we  shall,  in  general,  have  some  (n-l)-diraen8ional 
manifold,  on  which  the  value  of  T  is  prescribed.  Because  of  the  non-linearity 
of  the  eiconal  equation  (17)  there  may  be  several  solutions  with  these  prescribed 
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initial  values.  We  shall  think  of  each  solution  f  as  defined  only  on  one  side 
of  the  initial  manifold  (and  zero  on  the  other),  and  in  its  region  of  definition 
we  associate  with  f  a  corresponding  wave  (for  which  the  initial  values  of  the 
z    still  need  to  be  prescribed)  which  is  also  zero  outside  this  region.  Thus 
our  waveswill  in  general  be  discontinuous  across  their  initial  ma■^ifolds^ 

It  is,  by  now,  perhaps  obvious  what  must  b«  done  to  complete  the  con- 
struction of  V.  We  shall  attempt  to  repair  the  discontinuities  in  y  by  adding 
wavas  to  (10 ).  Thus  our  formal  solution  will  be 

a 
where  the  z    must  be  appropriately  specified.  Since  y  is  smooth  except  on  the 
manifold  P ,  we  take  P  to  be  the  initial  manifold  for  each  z   *  At  this  point 
vre  can  easily  state  the  outgoing  condition  referred  to  in  the  Introduction o  First 
we  note  that  with  every  solution  T,  and  therefore  with  every  wave,  is  associated 
an  (n-1)  parameter  family  of  rays,  A  wave  will  be  said  to  be  outgoing  (incoming) 
with  respect  to  the  initial  manifold  p  if  on  every  one  of  its  associated  rays, 
T  increases  (decreases)  monotonically  with  distance  from  p«  In  the  expansion 
of  (3U)  only  waves  z    which  are  outgoing  with  respect  to  p  are  to  be  included. 
The  requirement  that  v  shall  be  continuous  may  be  stated  simply  as 

(35)  pvjT  -  0    on  p, 

where  []^v~l'  denotes  the  jump  in  v  across  P  (i.e.,  the  difference  of  the  limiting 
values  as  p  is  approached  from  each  side)  with  some  given  orientationo  Inserting 
(3ti)  in  (35),  we  have 

(36)  r_Y2^   4  ^  e(a)  z^°^  -0   on  p . 

a 

n 
Here  e(a)  ■  *  1  and  it  is  determined  by  the  orientation  of  [  ] '  and  the  side  of 

P  on  which  z    is  zero.  Since  [  y] '  contains  no  exponential  factor  in  co  we 
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see  that  (36)  can  be  satisfied  only  if 

(37)  1^"^   -  0    on  r. 

Inserting  (6)  and  (13)  in  (36)   and  equating  the  coefficient  of  each  power  of  (i&)) 
to  7.ero,  we  obtain 

, (a)(0) 


(38)  J2^^^^^  -  0     on    r 


a 


(39)  E«(«)^^°^^'^  -  [^^^■^^g]'^,     on   Ti 


r  ■  i.^c^«*«     # 


Now  we  muet  determine  the  range  of  a,  that  is  we  must  determine  how  many 
outgoing  waves  appear  in  (3U).     In  solving  the  eiconal  equation  and  initial 
conditions  (37),  T        and  its  derivatives  must   satisfy  the  so-called    'strip 
manifold   condition'.     Geometrically  this  means  that  the  solution  surfao* 
T  ■  Y       (x)   in  Tx*spaoe  must  be  tangent  to    p ,   so  the  projection  into  x-space 
of  the  normal  to  the  solution  surface  must  be  parallel  to  the  normal  to    p ,     Thus 

(I4C)  T^<^^  -     piaW     i  V  -  l,...,n 

V  ^ 

where  Y-i,«»»,Y  are  components  of  the  unit  normal  to  p  in  x-space,  and  p(a)  is 

a  proportionality  factor.  It  is  easily  seen  from  the  eiconal  equation  that  a 

solution  surface  can  never  have  a  tangent  plane  parallel  to  the  plane  T  ■  0  (in 

Yx-space),  Hence  each  solution  f    is  increasing  (with  distance)  on  one  side  of 

P  and  decreasing  on  the  other.  Thus  each  T    provides  us  with  two  waves  one  on 

each  side  of  p .  One  is  outgoing  and  one  is  incoming.  We  use  only  the  outgoing 

cane.  The  scalar  factor  p(a)  is  determined  by  the  eiconal  equation  (17) J  thus 

(Ul)        det  G^*^^  -  det(l-  f~    T^^^A^)  =  det(l  -  p(a)E)  -  0,   on  P . 


6=1  ^v 


where 
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There  will  be   as  many  outg^oing  waves   as   there  are  distinct  solutions  p(a)   ■  •S7~T 
of  (bl),  hence  as  many  as  there  are  distinct  non-zero  eigenvalues  p(a)  of  E. 
(There  will  be  equally  many  incoming  waves,  namely  the  waves  associated  with  the 
continuation  of  each  f^°''  across    p.) 

Let  us  introduce  the  orthonormal  eigenvectors  B**  of  E  which  are  associated 
with  the  non-zero  eigenvalues  p(o)s 

(U3)  (I   -  p(j)E)B^     -     G^^^B^     -     0,         on   r. 

The  p(j)  and  G^*'     in  ()43)  are  not  necessarily  distinct  for  distinct  j,  i.e., 
p(j)  -  P(i)  and  g'^^   -  G^^^  if  p(3)  «  p(i).     Furthermore  each  non-mill  eigen- 
fector  B''  of  F  is  a  null  eigenvector  of  G^°     for  some  a,  and  the  totality  of 
null  eigenvectors  for  all  the  G         are  the  totality  of  non-null  eigenvectors 
B''  of  E,     This   statement  holds,  of  course  only  on   P  where  the  B     are  defined. 
The  matrix  E,  however,  may  W3ll  be  singular.     In  this  case  there  will  also  be 
null  eigenvectors  J   ,...,J*^  of  E.     In  order  to  comolete  the  construction  of  a 
continuous  formal  solution  v  it  is  necessary  to  impose  a  new  condition  on  the 
discontinuities  of  g(x)  which  we  call  the 
source  condition;  "^Ls]       "^»  v«l,...,p;       onp, 

(Here  we  have,  of  course,   an  inner-product  of  two  vectors.     If  E  is  non-singular 
the  condition  is  empty,  i.e.,  it  does  not  restrict  g.)     It  will  be  seen  in  the 
next  section  that  the  source  condition  has  a  very  interesting  interpretation  for 
Maxwell's  equations. 

We  must  now  proceed  to  select  initial  values  in  order  to  satisfy  (38) 
and  (39).     In  so  doing  we  will  make  constant  use  of  the  fact  that  the  set  of 
eigenvectors  B'^  may  be  decomposed  into  disjoint  subsets  each  of  which  is  precisely 
the  set  of  null  eigenvectors  of  one  of  the  G^°    .     Since  the  B''  are  independent  and 
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since  each  z^°'^  ^  has  an  expansion  in  terms  of  the  null  eigenvectors  of  G   ,  it 
follows  from  (38)  that  each  2^°'^  '  is  zero  on  P ,  Therefore  the  initial  values 
of  the  0^  ^°'     are  zero,  and  since  equation  (32)  is  homogeneous  for  r  »  0,  we 
conclude  that  z^°^^  ^  =  0  for  each  a  ,  We  now  consider  (39)  for  r  «  1.  Again 
each  z^"-^^   has  an  expansion  in  terms  of  the  null  eigenvectors  of  G   ,  for  by 
(2U)  0^^°^  -  0  because  z^°^'  ^5  0.  It  is  clear  now  that  ve  may  satisfy  (39) 
only  because  of  the  source  condition,  and  that  the  initial  values  on  P  of  each 
gVQA   ^^g  uniquely  detemaned.  Now  suppose  we  have  determined  the  z^^'^^     and 
have  satisfied  (39)  for  p  ■  0,l,,..,r-l.  We  recall  that  z^'^^^'*^  has  an  expansion 
in  terms  of  all  the  eigenvectors  of  G   .  By  (2U)  the  coefficients  of  the  ncn- 
null  eigenvectors  are  determined  in  terras  of  2^°'^    ,     It  remain?  to  determine 
the  coefficients  of  the  null  eigenvectors  from  (39).  This  can  be  done,  and  (39) 
will  be  satisfied  if  and  only  if 


ihh) 


^      f    s.v   (a)(r)   .vpy(r-l)"ir      p      . 

2_  e(a)J  z      •  J   c^  '   'g  '  ,  on  Pj  v  -  l,...,p. 


We  will  now  show  that  (UU)  ie  a  consequence  of  the  foregoing  induction 
hypothesis.  First  we  consider  the  directional  derivative 

for  any  v  -  l,.,.,p.  It  is  easily  seen  that  the  direction  of  —  is  tangential 
to  p«  Hence  for  any  vector  u  which  is  smooth  on  both  sides  of  p. 


(U6) 


Hence 


d 
dr 


wr.[3^q^. 


(1.7)    jV [u]  r.  I-  [,]  r,  jVgj.^^ r,  ^^  „-j  r.  ^^^^  r.  pv^^j  r .  /j-^  q  r 


It  follows  that 


i.e.,  we  could  have  started  the  expansion  (13)  with  r  «  1 
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(W) 


jV[5^^-^^]r.4/(-^)g]r. 


We  observ©  that  each  J     is  a  non-null  eigenvector  for  every  G       ,     Thus  from  (2U), 

But 

(50)  jV^^j"     .     j"  (I-  p(a)E)j''     -     1. 

Hence 

(51)  jV°)(^)     .    jVz^^^^"*"^^  . 

By  the  induction  hypothesis  we  may  assume  the  validity  of  (39)  with  r  replaced  by 
r-l,  i.e.. 

Now  using  (51),  (52),  and  (b8),  (in  that  order)  we  have 

(^3)        Ee(a)jV°^^^^  -  jV  E«(a)^^'^^^"'^-  j"<=2fp!^""^^gF 
a  a  L     J 

Thus  we  have  verified  (Uii)  and  the  induction  argxunent  is  complete. 

It  follows  that  the  initial  values  for  z      will  be  uniquely  determined 
for  all  a  and  r  and  (39)  will  be  satisfied  for  all  r.  We  may  then  solve  the 
transport  equations  recursively  for  r  -  0,1,,,,  along  each  ray,  starting  at  p , 
and  determine  the  z      ,  With  this,  the  formal  construction  is  complete  except 
for  the  following  reservation:  Examination  of  the  transport  equations  reveals 
that  the  coefficients  of  the  linear  differential  equation  (32)  depend  on  the 


-  hi  - 

eigenvectors  of  the  matrix  G  ,  hence  on  the  function  T^°  ,  By  our  assumptions 
about  p,  1  will  be  smooth  in  a  neighborhood  of  p,  but  as  we  proceed  along  a 
ray  from  p  we  may  reach  a  point  where  T  or  one  of  its  derivatives  becomes 
singular.  We  cannot,  by  the  present  construction,  continue  the  solutions  of  the 
transport  equations  beyond  such  a  point.  Points  of  this  type  are  said  to  lie  on 
a  caustic  •  For  convenience  we  state  a  condition  which  guarantees  that  the  con- 
struction can  be  completed: 

caustic  condition;     The  problem  is  such  that  the  phase  functions  T       ,   associated 
with  each  of  the  waves  in  the   cwistruction,   are  free  of  caustics. 

Having  completed  the  formal  construction  for  the  probUero  as  formulated 
in  the  Introduction,  we  shall  now  show  how  to  extend  the  construction  to  problems 
involving  boundaries. 

If  p(x)   •  0  is  the  equation  of  any  smooth  surface  B  in  the  x-space,  the 
transformation  z     ■  a  (x);  v  =■  l,...,n;  with  a     •  p  maps  B  into  the  hyperplane 
z     ■  0,  while  preserving  the   syirmetri  c  hyporbolic  nature  of  the  operator  L,  and 
therefore  the  properties  of  the  reduced  operator  N.     Thus,  without  loss  of 
generality  we  may  take  our  boundary  B  to  be  the  plane  x     -  0.     On  this  boxindary 
we  require  our  solution  v  (hence  our  for:nal  solution  v)  to  satisfy  the   following 
type  of  homogeneous  boundary  condition. 

The  vector  v  is  an  element  of  a  k-dimensional  Euclidean  vector  space  V, 
Consider  a  subspace  T  with  the  following  properties: 

It,     The  matrix  a"  is  non-positive  over  T,   i.e.,  uA  u  <  0   for  every  u  in  T. 

2,     T  is  maximal ,  i.e.,  the  dimension  of  T  is  as  large  as  the  dimension  of 
any  subspace  having  property  1. 

Our  boundary  condition  may  now  be  expressed  simply  as 

(5li)  v{x)      in     T,       for     x^     »     0, 


-a- 
For  a  discussion  of  the  formal  construction  for  certain  problems  involving  caustics. 


see 


M. 
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If  T"^  is  the  orthogonal  complement  of  T,  and  T"^  is  generated  by  the 
vectors  b  ,«..,b  ;  the  condition  may  also  be  stated  as 

(55)  vb     ■  Oj  3  ■  l,,..,in}  for     x     ■  0, 

y 

The  vectors  b     are  assumed  to  be  smooth  fxmctions  of  x^.,.*,x     ., 

1'   '  n-1 

For  problems  with  boundaries,  we  have  assumed  that  the  source  region 
does  not  intersect  the  boundary,  i.e.,  that  g(x)  «  0  in  a  neighborhood  of  the 
boundary.  Thus  the  term  y  in  the  expansion  (20)  of  v  automatically  satisfies 
the  boundary  condition.  Suppose  now  that  a  particular  wave  z    (outgoing  from 
P)  is  'incident'  on  the  boundary,  B  (i.e.,  any  of  the  associated  rays  intersect 
the  boundary).  We  observe  first  that  z    will  be  incoming  with  respect  to  B. 
Our  procedure  will  be  to  add  appropriate  reflected  waves  to  the  expansion  v 
which  will  be  outgoing  with  respect  to  B,  and  which,  together  with  z^   will 
satisfy  the  boundary  condition.  We  emphasize  that  the  condition  will  be 
satisfied  individually  for  each  z    together  with  its  associated  reflected 
waves.  We  shall  see  that  the  boundary  condition  uniquely  determines  the  re- 
flected waves. 

The  incoming  wave  z    will  have  some  given  phase  (i.e.,  the  value  of 

T   )  on  B.  It  is  clear  that  the  associated  reflected  waves  must  have  the 

same  phase  on  B.  In  geometrical  language,  if  the  phase  surface  Y  «  T   (x) 

in  Tx-space  intersects  the  boundary  B  (i.e.,  the  manifold  x  «  0  in  Yx -space) 

in  the  (n-l)-diraensional  manifold  ^,  then  the  corresponding  phase  surfaces 

of  the  reflected  wa'vos  must  also  intersect  B  in  /\,    We  shall  sometimes  denote 

by  x  the  coordinate  Y  in  Yx-spac«. 

Let  5  ,...,6  ,,0  denote  components  of  the  unit  normal  to  /\  in  the 
o'     n-l 

plane  x  -  0,  Let  Y  •»  Y^  (x)  be  the  equation  of  a  phase  surface  C  '^  whose 
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intersection  with  B  is  /\,  Then  (by  the  strip  manifold  condition),  since  the 
ccniponents  of  the  normal  to  C  '^  are  (1,  -T^  ,,,.,-Y  "^  ),  we  have  1  ■  p6  } 
-T^P^  -  p5  ,  V  =  l,...,n-l:  and  T^^^  i 


■  p5  ,  V  =  l,..,,n-lj  and  y^*^'  is  a  solution  of 
*v  n 


(56) 


Her* 


det  G^P^  -  detd  -  r  T^P^A^)  =  det(-T^P)A"  ♦  pF)  =  C,  on  Z^  . 


— ^      X  'X 

v*l   V  n 


n-1 

(57)  F  -  ^   6  A^. 

v»0 

Equation  (56)  will,  in  general,  have   complex  roots  f  ^   ,     For  each  distinct 

n 
real  root  there  will  be  a  corresponding  phase   surface  passing  through  /^  • 

We  now  impose   a  new  condition  on  the  problem,  the 
shadow  condition!     F  shall  be  non-singular  throughout  B, 

In  Part  I,   Section  6     it  is  pointed  out  that  at  the  point  of  rainimujn 
T       ,  F  ■  I,     Therefore  it  follows  from  the  above  condition  that  F  is  positive 
definite  on  B«     It  is  also  pointed  out  that   the  condition  guarantees  that  no 
rays  of  either  the  incoming  or  reflected  waves  can  be  tangent  to  the  boundary. 
Thus  we   can  have  no  shadows  in  the  problem. 

Since  F  is  positive  definite  we  may  set  F  ■  JJ'  where  J  is  a  non-singular 
matrix  hence  has  an  inverse,  S,     Let    au    =  SA  S'.     Then     kj     is  symmetric  and 
we  may  write 

(58)  (;^  -  Op)   Q^P^   -  Oj  p  =  l,...,k 

Here  the  Q^^'  are  a  complete  set  of  orthogonal  eigenvectors  of  )0  ,  which  we 
may  normalize  by  the  condition 

(59)  Q^P^SS'   Q^P^  -  Ij       p  -  l,...,k} 

The  eigenvalues  a  are  real.  Set  B^^'   -  S'q'P'j  p  =  l,,.,,k.  It  follows  from 

The  discussion  of  that  section  closely  parallels  this  one  if  we  identify  T  with 
t,  'phase  surface'  with  'characteristic  siirface',  etc. 
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(59)  that  the  B^^  are  unit  vectors.  They  are  clearly  linearly  independent. 

From  (58),  J(Sa"s'  -  a  )J'B^P^  -  0,  or 

P 

(60)  (a""  -  a  F)B^P^  -  Oj      p-l,...,k. 


For  each  non-zero  a  ,   set  1^^'m  p/a  ;  then 

P  3t^  P 

(61)  G^P^B^P^   -   (-T^P\"  ♦  pF)B^P^   -  0;        U  a^O. 

n  ** 

Equation  (61)  expresses  the  fact  that  for  each  distinct  non-zero  a     there  is 

P 

a  distinct  phase  surface,  C  P  ,  passing  through  Zl,  and  B  P  is  a  null  eigen- 
vector of  the  associated  matrix  G^P\  Since  we  may  have  multiple  roots  o  , 
we  may  have  several  distinct  B^P"^  associated  with  tte  same  matrix  G^P', 

We  have  up  until  now  not  distinguished  between  waves  2  P  which  are 
incoming  and  outgoing  with  respect  to  B  .  For  convenience  we  simply  assign 
the  value  zero  to  the  coefficients  of  all  incoming  waves  except  the  one  z^*^' 
which  "  arises"  at  p.  We  will  then  attempt  to  determine  the  initial  values 
for  the  outgoing  waves  from  the  boundary  condition: 

(62)  b^  Yl    *    "  °>  °^^»  ^  ■  !*•••»">• 

P 

Comparing  powers  of  (ico)  in  the  fxpansions  of  z^P'^  it  follows  at  once  that 

(63)  b^  ^z(P)(^)  -0}  onB}     ""   ^'•••'"' 

p  r  ■  0,1,2,... 

We  recall  that  each  z^'P''^  '   has  an  expansion  in  terms  of  null  eigenvectors  of 
G^P  ,  But  this  set  of  vectors  is  just  a  subset  of  the  set  of  vectors  B  "^  . 
Thus 

"5 ' 

We  are,  of  course,  dealing  only  with  waves  on  one  side  of  B,  i.e.,  within  the 

region  of  the  problem,  which  we  take  to  be  x  >  0. 
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(The  O"  may  be  identified  with  the  6"  in  (18),  but  ve  have  renumbered  them  and 
p  m 

(6U)  contains  the  coefficients  for  several  different  waves.)  Inserting  (6U)  tn 
(63)  for  r  -  0, 

(65;)  Tb^B^PV^-Oj     v-l,...,in. 

P^  ^ 


Of  course,  some  of  the  C~^   are  already  determined,  namely  those  associated  with 

o  ^^^^ 
the  incoming  waves.  It  is  easy  to  see  that  6" 


will  not 


value  at  a  point  on  B  if  Y  U  ^r    with  x  ,   i,e.»  if 

'^  jincreasesj  n'  * 


have   a  preassigned 


dx         dx  /    \       /    \ 

i    -     'i  -     "  -  b(p^a"b(p^    i. 

n       ds         dT 


<  0 
>  0 


(The  possibility  x^  -  0  is  precluded  by  the  shadow  condition.)     Let 
b(p'a"b(p)  b,   |>  °  III  P  :  l-;.;=^J  .    „e  rewlte  (65)  pUoine  known 
quantities  (which  we  denote  simply  by  W  )  on  the  right  side: 

(66)  IIb^B^PV°-W;  V  -  l,...,in. 

p«l  ^ 

In  Part  I,  Section  6  it  is  proved  that  the  coefficient  matrix  (b'^B^P')  of  (66) 
is  square   (i.e.,  s  »  m)  and  non-singular.     Hence  (66)  can  be  solved  uniquely 
for  <r°,...,(r°j  and  these  quantities  provide  the  Initial  values  for  the 
transport  equations  of  order  zero   (r  -  0)  for  the  reflected  waves. 

To  obtain  the  initial  conditions  for  the  higher  transport  equations, 
we  use  (63)   for  r  »  1,2,..,      .    We  recall  that  each  z^^^^     has  an  expansion 
in  terms  of  the  eigenvectors  of  G^^^  the  coefficients  of  the  non-null  eigen- 
vectors being  known  from  (2U)  in  terms  of  z'^'^^'     ,     Thus  we  may  write 

(67)  y:  Z^P^^^^    .   E   ^^B^P^  .W2  . 

p  P-1 


-  52  - 

From  (63) 

k  -   /  ^ 

(68)  T"  b^B^P^  <r^  =  W,,    V  .  l,...,in. 

As  before,  those  C     associated  with  incoming  waves  are  known.  Thus 

P 

(69)  F  b^B^P^  tf'J'  «  W,,,     v-l,...,ni.- 


& 


Since  (66)  and  (69)  hsve  the  same  coefficient  matrix  it  is  clear  that  the 
latter  also  has  a  unioue  solution.  Thus  the  initial  values  for  the  transport 
equations  of  all  orders  (for  the  reflected  waves)  are  uniquely  determined. 

Assuming,  now,  that  the  caustic  condition  is  satisfied,  we  may  solve 
the  transport  equations  and  cojr.pletely  determine  all  of  the  reflected  waves 
(for  all  of  the  incident  waves  2^°').  This  completes  the  construction  of  the 
formal  solution  for  problems  with  boundaries  «  If  we  now  use  the  index  a  to 
label  all  our  waves,  including  reflected  ones, the  formal  expansion  may  be 
written  in  the  form 


(70)  V- 


^      f.    Tr+l    ^  ^   /.  sr 

r«c  (ico)      a  r-1   (iw) 


2c     Maxwell's  equations 

In  the  present  section  we  will  demonstrate   the  applicability  of  the 
results  of  this  paper  to  Maxwell's  equations.    We  begin  by  stating  the  equations 
in  their  usual  form: 

(71)  VxH  -1|^  (eE)  -  i^Ij  I  -  <rE  +  ^F(x,t) 

7 — 

For  problems  with  multiple  boundaries  and  multiple  reflections  the   construction 

is  an  obvious  extension  of  the  foregoing. 
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(72)  VxE  +  i|f.(MH)   .  0 

(73)  V   •  (eE)  =  Utip 
(7U)     V  •  (pH)  -  0. 

Here  the  independent  variables  are  t  *  x     and  x  «  (x^,x_,x-);  H,  E, 
I,  F  are  vectors  with  three  ccmponente.  The  scalars  e,  ix,  a"   will  be  taken  to 
be  sBooth  functions  of  xj  e  and  ^i  are  assumed  to  be  everywhere  positive. 
(Physically  we  are  restricting  consideration  to  isotropic  media  in  which 
the  electromagnetic  properties  do  not  change  with  time,  and  vary  smoothly 
with  position.)  F  represents  impressed  currents,  i.e.,  sources. 

We  shall  show  that  equations  (71)  and  (72)  form  a  symmetric  hyperbolic 
system.  As  such,  the  usual  existence  and  uniqueness  theorems  for  initial  value 
problems  (see  Part  I)  api:)ly.  In  view  of  the  completeness  of  this  theory 
equations  (73),  (7h)  appear  to  be  somehow  superfluous.  The  confusion  is  easily 
resolved,  however,  if  we  interpret  them  as  additional  conditions  on  the  initial 
data  of  the  problem.  From  this  point  of  view,  we  begin  with  (71)  and  (72); 
then  we  define  the  fimction  p(x,t)  by  equation  (73),  and  impose  the  following 
condition  on  the  initial  data: 

(75)  V   •  (nii)  «  0   at   t  -  0. 

It  follows  easily  from  (72)  that 

(76)  |-  (V  .  tjH)  -  0, 

and  ilk)   is  an  immediate  consequence  of  (75)  and  (76).  From  (71)  and  (73)  we 
obtain  the  'equation  of  continuity', 

(77)  7  •  I  +  p^  -  0. 
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We  now  set  E  -  (0^^ ,02*12^3),  H  -  i<^i^,^^»0^) ,  F  -  (f^,f2,f3). 


Gf. 


f     - 


0 
0 
0 


Then  (71)  and  (7?)  may  be  written  as 


(78) 


3 


C°  0,   +     E       C\     ♦     ^  0  *  f 


v-1 


where  the  matrices  C°  and  /^    are  given  by: 


(79) 


c°  -  5^  3- 


1 


/I  0  0  0  0 
0  /|  0  0  0  0 
0  0  )^  0  0  0 
0  0  0  /|  0  0 
0  0  0  0  /I  0 
0     0     0     0    0    /I 


,^.!!i£- 


V 


1^1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

and  the  matrices  C     may  be  conveniently  written  together  in  the  form: 


0  0  0  0 
0  0  0  -z 
0       0 


^  ■-  n  v' 


v»l 


Z3     -Z2 


I 


0       Zg     -Zj^     0 


0      -z,     Zg     0       0         0 

z       0     -z^     0       0         0 

-Zg     z^     0       0       0         0 
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-1. 


To  transform  (78)  into  our  standard  form,  we  set  0--6'     u,     Ifwe  insert 
this  expression  in  (78)   and  multiply  on  the  left  by     <[       we  obtain 


(80) 

where 

(81) 

(82) 

(83) 


(8U) 


\    * 


V"    a\      +  Bu 

(■ — ^_    X 


v-1 


^— ,  X  ^— -  X 

V=l  V  V=l  V 


n 


-  ^-If. 


It  is  clear  from  either  (73)  or  (80)  that  we  have  a  symmetric  hyperbolic 
system.  Furthermore  if  we  define  the  matrix 


3 


v»l 


a''  ^  3'^m 


-1  n/  ^  -1 


■  e 


^; 


_c 


one  may  verify  that  the  eigenvalues  of  P  are 
(85)         0,  0,  ed,  ed,  -ed,  -ed  ; 


/? 


^ 5 ? 

+  Z2  +  z-   • 


It  follows  at  once  that  the  multiplicity  condition  is  satisfied. 

For  convenience,  we  also  list  a  set  of  orthonormal  eigenvectors  of  P 
in  the  same  order  as  their  corresponding  eigenvalues  appear  in  (85): 
(86) 


d  /? 


R 


d  /? 


-z. 


-z. 
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(86  cont'd.) 


i>4u|r^ 


- -2^23  - 

■'2=^3 

2     2 

dz^ 

-dz^ 

-  0  ^ 

1/2(2^*2^) 


dZp 
0  ^ 


2^23 
'2^3 


f   2     ?x 
-(z^+z^) 


m-fyT^) 


2^2^ 
2^2^ 

,'    2      2. 

2,d 

-z^d 


iyfczJTzf) 


d22 
-dz^ 

0 


-2^23 

-^2^3 

2     2 

V^2 


For  Maxvrell's  equations,  we  shall  take  our  boundaries  to  be    'perfect 
conductors',  i.e.,  we  shall   apply  the  boundary  condition 


(87) 


■"tan 


0 


where  F    denotes  the  components  of  E  tangential  to  the  boundary.  If  we 
denote  the  boundary  by  y^(x)  »  0,  the  transformation  y  »  y  (x);  v  »  1,2,3  leads 

to  an  equation  in  which  the  coefficient  matrix  A_  multiplying  u   is  simply  the 

ni  ^^3  ^^ 

matrix  P  -  e  yy  ,  provided  we  set  z  =  ~-  •  Since  n  »  (2,  ,Zp,z,)  now  denotes 

a  vector  normal  to  the  boundary,  the  condition  may  be  expressed  as 


(88) 


u.  ■  ca.; 


J  =  1,2,3 


where  u.  denotes  components  of  u.  It  is  clear  that  all  -"sctors  n  satisfying 
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(88)  lie  in  a  space  T  of  dinension  k»     Furthermore,  for  u  in  T  it  is  easily 
verified  that  ]i  A  u  «  0,     Thus  ^     is  non-positive  over  T.     Since  i    -  P  has 
precisely  h  non-positive  eigenvalues,  it  is  easily  shown  (set*  Part  I,  Appendix 
2)  that  T  is  maximal.     Thus  we  have  shown  that  the  boundary  condition  (8?)  is 
of  the  form  we  have  applied  in  Section  1. 

We  examine,  now,  the    'source  condition',  vhich  was  introduced  in 
Section  1,     If  y  "  (Yi»Yo»Yo)  denotes  a  unit  vector  normal  to  the  surface  of 

discontinuity    p   i"  ^^   source  region,   and  J         are  null  eigenvectors  of  the 

^  V 

matrix  E  -  2^    Y  A  ,  this  condition  is 

V"l      ^ 

(89)  «J^Cg  ]  -0       (for  all   J^). 
The  function  g(x)  is  determined  by 

(90)  ge-i^  .  >j    .  -    .^  -1  f  -  -    /I     f . 

We  may  identify  the  matrix  E  with  the  matrix  P  by  setting  y    •  a   j   v  -  1,2,3. 
Thus  the   JT  are  given  by  R     and  R     of  (86).     Since  e  is  a  continuous  function, 
we  may  rewrite   (89)  as 

R^^Qf  ]     -    Oj  11  -  1,2; 

or  in  terms  of  3  component  vectors,        ,  " 

(91)  yQf]  -  Q^]  -  0. 

Thus,  in  physical  terras  the   source  condition  may  be   stated  as  follows:     Across 
any  surface,  the  normal  component  of  the  impressed  current  must  be  continuous. 
In  Section  6  we  shall  introduce  the  so-called  initial  data  condition  : 

(92)  J^[i(0,x)1   -  0  (for  all  J^"). 
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For  Maxwell's  equations  this  becomes 
(93)  R^pO,x)1   -     0;         ti  -  1,2    . 

Equation   (93)  is  equivalent  to  the  two  equations 

(9I4)  [yH]     -     0       at     t  -  Oj 

(95)  QyE  ]     -     0       at     t  -  0. 

Equation  (9U)  is  a  direct  consequence  of  (7U),     Equation  (95)  imist  be 
interpreted  as  an  additional  condition  on  the  initial  data  of  the  problem.     In 
physical  terras,   (95)    states  that  there   shall  be  no    'surface   charges'  present, 

SiMmarizing  the  results  of  this  section  we  may  say  that 

1)  Maxwell's  equations  are  sym^nstric  hyperbolic  and  satisfy  the  multiplicity 
condition . 

2)  The  boundary  condition  for  perfectly  conducting  boundaries  is  a  special 
case  of  the  type  we  considered  in  Section  1, 

3)  The  soiirce  condition  is  satisfied  for  impressed   currents  which  satisfy 
(91). 

U)     The  initial  data  condition  is  a  consequence  of  the  divergence  equations 
(73),   (7b)  provided  we  exclude  surface   charges. 

It  follows  from  1),  2)  and  3)  above  that  our  formal  expansion,  and  the 
entire  contents  of  this  paper  are  applicable  to  Maxwell's  equations.  The  signi- 
ficance of  U)  is  discussed  in  Section  6. 
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3.   The  Duhamel  integral  theorem 

In  the  introduction  to  Part  II,  we  stated  that  the  results  of  Part 
I  play  a  key  role  in  the  proof  of  the  asymptotic  nature  of  our  formal 
solution.  In  Part  II  w©  have  been  concerned  with  an  inhomogeneous  hyperbolic 
equation  (2)  (and  its  related  reduced  equation  (U))  whereas  in  Part  I  we  were 
concerned  with  the  homogeneous  equation  Lu  »  0 ,  The  connection  between  these 
two  equations  is  established  by  the  !>ihamel  integral  theorem.  Before  proving 
this  theorem,  however,  we  must  discuss  the  notion  of  'veak  solutions'. 

As  in  Part  I,  we  consider  a  space  ^q   of  smooth  testing  functions 
_(\  (t,x)  with  compact  support.  More  precisely  ^r\      '^'^^'^   ^®  '^^®  collection 
of  all  smooth  k  component  vector  functions  J^  each  of  which  vanishes  identi- 
cally outside  some  compact  region  P.  which  depends  upon  XL  and  lies  in  the 


half  space  t  >  0,  (For  problems  with  boundaries  we  assume  that  R  does  not 
intersect  the  boundary,)  For  any  u  and  f  in  C  ,  the  set  of  functions  with 
continuous  first  derivatives,  and  any  J\  in     W/-\  we  have  the  identities 

(96)  I  il.  Lu  dxdt  »     u  M  XL  dxdt}     _fL  f^<ixdt  -  "  j  f  i\+  dxdt. 
Here  the  operator  M,  the  adjoint  of  L  is  defined  by 

(97)  Mil  --X:   (A^-TL)^  *  B'  XL  . 

v«0        '^v 

B'  denotes  the  transpose  of  the  matrix  B.  The  integration  in  (96)  extends  over 
the  half  space  t  >  0.  The  identities  are  obtained  simply  by  integration  by 
parts,  the  boundary  terms  vanishing  because  X\.  has  compact  support.  If  now 
we  have  Lu  »  f . ,  it  follows  trivially  that 
//  / 


(98) 


uM  XL  dxdt  + 


f  XL^  dxdt  -  0,   for  all  XI  in  ^^» 


Conversely  if  u  and  f  are  in  C  in  the  neighborhood  of  some  point  in  the  half 
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space  t  >  0  and  satisfy   (98)  it  follows  that  Lu  ■  f .  in  that  neighborhood. 
VJe  may  however  consider  more  general  classes  of  functions  than  C  .  Thus  if 
u  is  any  integrable  function  which  satisfies  (98)  we  shall  say  that  u  is  a 
weak  solution  of  the  equation  Lu  =  f.  .  If,  in  addition,  u(0,x)  ■  g(x)  (and 
u  is  in  T  on  X  ■  O)  we  say  that  u  is  a  weak  solution  (I)  of  the  problem 

(99)      Lu  -  f^j   u(0,x)  -  g(x),     (u  in  T  on  x^  =  0). 


VJe  shall  also  need  a  second  definition  of  weak  solutions.  For  this 
purpo5e  we  consider  a  second  space  Jo-m    of  smooth  testing  functions  J  ,  w^iich 
have  the  same  properties  as  the  functions  J~L  except  that  we  now  require  only 
that  the  function  J  vanishes  identically  outside  some  region  R  whose  inter- 
section with  the  half  space  t  >  0  is  compact  (see  Fig.  3)»  (For  problems 
mth  boundaries  we  again  assume  that  R  does  not  intersect  the  boundary.)  It 
is  clear  that  S ^^  is  a  subset  of  S^,    We  now  consider  the  analogue  of  (98) 


(100) 


u  M  I  dxdt  +   M"  J  t,  ^^^^  '   °      ^°^  ^^"^  ^  ^"  ^IS 


As  before  if  u  and  f  are  in  C  in  a  neighborhood  of  any  point  in  the  half  space 
t  >  0,  then  Lu  »  f .  in  that  neighborhood.  If  now,  u  and  f  are  in  C  in  a  neigh- 
borhood of  a  point  (0,x)  on  the  initial  plane  t  =  0,  we  may  integrate  by  parts 
to  obtain 

■5 "" 

This  is  the  boundary  condition  defined  in  Section  1.  In  the  fol] owing  sections 

we  shall  frequently  treat  initial  value  problems  and  initial  boundary  value 

problems  simultaneously.  Parentheses  are  used  to  denote  the  additional 

conditions  or  results  associated  with  the  problem  with  boundaries. 
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(101)       0  -    I  I  J  (Lu-f  )cbcdt  + 


I       J(0,x)Ju(0,x)-f(0,x)ldx 


f(0,x)^u(0,x)   -  f(0,x)V    dx. 
't-0  [  J 

It  follows  that 

(102)  u(0,x)     -     f(0,x). 

If  u  and  f  are  any  integrable  functions  which  satisfy  (100)  and  f. 
exists  (and  u  is  in  T  on  x  =  0)  we  say  that  u  is  a  weak  solution  (II )  of  the 

problem 

(103)  Lu  -  f. }    u(0,x)  -  f(0,x);     (u  in  T  on  x^  -  0)  , 

In  Part  I  we  introduced  the  class  }\^   of  piece-wise  smooth  functions: 
For  each  vector  function  u  in  /{  t   defined  in  the  region  of  interest  ^       , 
there  are  a  finite  number  of  smooth  hypersurfaces  which  divide  /<;  into  subregions 

fQ      such  that  the  first  derivatives  of  u  can  be  made  continuous  in  the  closure 
of  each   ^  by  appropriate  definition  of  the  derivative  on  the  boundary. 

Our  two  definitions  of  weak  solutions  are  related  by  the  following 

Lemma!  If  u  and  f  are  in  ^,  then  u  is  a  weak  solution  (l)  of 
the  problem  (103)  if  and  only  if  it  is  a  weak  solution  (II)  of 
the  same  problem. 

Proof I 

1«  Let  u  in  /(  he  a   weak  solution  (II ),  Then,  as  before,  u(0,x)  -  f(0,x}. 
Since  S    is  a  subset  of  S*  ,  (98)  holds.  Hence  u  is  a  weak  solution  (I)  . 

2,  Let  u  in  ^  be  a  weak  solution  (l)  ,  and  let  C  denote  the  hypersurfaces 
subdividing  -C/  (for  both  functions,  u  and  f).  Then  for  any  J  in   »/» 

By  setting  f  •    q(x,t-)dr'  +  i  (x) ,   we  may  define  a  weak  solution  (7?)  of  the 

/  o 
initial  value  problem  Lu  =  q(x,t)}   u(0,x)  ■  /(x), 

Hit 

^   is  either  the  half  space  t  >  0,  or  in  problems  with  boundaries,  that  portion 
of  the  half  space  in  which  the  solution  is  defined. 
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(loU) 


uM  ^  dxdt 


f  f^dxdt  -     j  J  Jlu-f^l  dxdt  +    j  |(x,0)Ju(0,x)-f(0,x)  >■ 


6x 


where 


^  ,u,f 


denotes  an  integral  over  C     involving  the  jumps  in  u  and  f.     As  we  have  seen 
before,  Lu  ■  f     except  on  the  C    .     Also  u(0,x)  »  f(0,x),   hence  the  first  two 
terms  on  the  right  side  vanish.     Thus 


(105) 


u  M  I  dxdt  +  f  J     dxdt     -     J^   \    "<   J  »">f 


For  any     Pl  in       J  r\ 


(106) 


u  M  Si.  dxdt     + 


XV^dxdt  -  ^        J   il  ,  u,   fl  . 


/ 


It  is  clear  that  we  may  choose  an  S\.  which  is  equal  to  J  on  the  hypersurfaces 
C    .     Since  the  left  side  of  (106)  is  zero,  u  must  satisfy  (100),     Hence  u  is 
a  weak  solution  (II ). 

We  come  now  to  the  Duhanel  integral  theorem.     The  theorem  was  deduced 
by  R,  K,  Luneburp,  for  Maxwell's  equations,  using  Laplace  transforms,     M,  Kline 
proved  it  for  piecewise  smooth  solutions  of  Maxv^ell's  equations  and  the  second 
order  linear  hyperbolic  equation  by  a  process  of  direct  substitution  which  is 
quite  laborious.     The  following  version  is  both  more  general  and  more  direct; 

Theorem:     Let  u  be  a  weak  solution  (ll)  of  the  problem 


(107) 


Lu 


Oj         u(0,x)  -  g(x)j 


(u  in  T  on  X     »  0)o 
n 
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(108) 


Let  h(t)  be  a  smooth  scalar  function  of  t  and  set 
t 


W(t,x)   -  1^      I     u(s,x)h(t-s)ds, 


Then  W  is  a  weak  solution  (II )  of  the  problem, 


(109)         LW  -  g(x)h'(t);   W(0,x)  -  g(x)h(0)}    (W  in  T  on  x^  -  0) . 
Proof;  If  we  consider  the  boundary  condition  as  in  ($5)  it  is 
easy  to  see  that  W  satisfies  the  condition  if  u  does.  Now  by  definition 


(110) 


uM  J  dxdt  +         g(x)  J^   dxdt  -  0  for  all  J  in     J  *  , 


Usinp  integration  by  parts  we  have 

,00  ,00  t 


/ 


(111)  'M  J  dt 


u(s,x)h(t-s)d8|  M  f  (t,x)dt 


Bt 

""'o 

00  t 

I      M  u(s,x)h(t-s)ds"|    |M  I  (t,x)J^ 

0  0  ~ 
00                     00 

1  u(s,x)    [   [m  I  (t,x)Th(t-s)dtds 

1?  s 

00  00 


dt 


-   j  u(t,x)    j      t  J  (s,x)T    h(s-t)dsdt. 


let 
(112) 


00  00 

^  (t,x)  -  -        ^^(s,x)h(s-t)ds  -  -        J^(7,+t,x)h(z)dz 


It  18  easily  seen  that  ^  is  in      </  7K  •     Also 


I 


00 


.00 


(113) 


M-f    -  -   I      Im  J  (z+t,x)T    h(z)dz  -  -   I    |m  J(s,x)j^h(s-t)ds. 
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(Compare  the  last  terms  of  (113)  and  (111),)     Now 


(IIU) 


00  /^  /°°  00  00 

f^(t,x)dt  «  -  J^^(z+t,x)h(z)d2dt  «  -      h(t)|  J^Jt+z,x)dzdt 


o     o 

,00 


o  o 


h(t)  f^(t,x)dt. 


From  (111),   (113),   (110),  (llh). 


(115)      I  WM  J  dxdt  -    u  M  -J  dxdt  -  -   g(x)   f^  dtdx 

g(x)h(t)  J^  dtdx, 

)J 

and  the  proof  is  complete. 

In  what  follows  we  shall  deal  only  with  functions  u,  w,  f  which  are 
in   ^,     By  the  lemma,  we  may  without  arobigtdty  refer  simply  to    'weak  solutions', 
Setting  h(t)  ■  -e~       /ico  we  obtain  the  following 

Corollary t     Let  u  in   ^  be  a  weak  solution  of  the  problem 


(116) 


Lu  »  Oj       u(0,x)  ■  g(x);  (u  in  T  on  x^  •  0) 


Set 


-icot 


(117)  w(t,x)  -  ^   [^^V-     J     "(«>^)«^'^  ^^ 


u(t,x)         -icot     (      /       \   i'«^j- 
^ .*   '  +  e  u(s,x)e       ds. 


ICO 


Then  w  is  in    ^  and  is  a  weak  solution  of  the  probL 


em 


(118) 


-  g(x)   . 


Lw  -  g(x)e"^^j     w(0,x)   -  -  -Sgi  ;  (w  in  T  on  x^  -  0). 
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In  the  remainder  of  this  paper,  u  and  w  will  denote  the  functions 
defined  by  (116)  and  (117),  where  g(x)  is  of  type  g.  Theorem  II -B  of  Part  I 
asserts  that  under  appropriate  conditions,  which  will  be  satisfied,  the  solution 
u  of  (116)  exists  and  is  unique  in  ->v  .  We  have  shown  above  that  (117)  is  a 
solution  of  (118).  The  uniqueness  of  the  solution  of  (118)  follows  easily 
from  the  uniqueness  assertion  of  Theorem  II -B,  In  order  to  avoid  confusion 
it  is  useful  to  remember  that  the  function  g(x)  plays  three  roles;  It  is  the 
initial  data  for  u,  a  factor  in  the  inhcmogeneous  term  of  the  equation  for  w, 
and  is  the  inhcmogeneous  term  in  the  reduced  equation,  Nv  ■  g. 

Ij.   The  steady  state  function 

We  begin  by  imposing  a  new  condition  on  the  operator  L  which  guarantees 
that  the  function  u(t,x)  decays  to  zero  as  t  — >ot): 
decay  condition;  If  u  is  the  solution  of  the  problem 

(119)         Lu  -  0;     u(0,x)  -  g(x)}       (u  in  T  on  x^  •  O) 

whpre  g(x)  is  any  function  of  type  g,  there  is  a  constant  oo    >  0  such  that  for 


every  x  in     /(J  ,  and  for  co  >  go    the  integral 


—    o 

00 


/ .      V   ioat . . 
u(t,x)e       dt 


exists  and  is  bounded  independently  of  oj,     . 

Since  for  any  n,  u       is  the   solution  of  a  problem  of  the  type   (119)  it 
t" 

follows  immediately  ftrom  the  above   condition  that 

In  the  work  of  Luneburg  and  Kline,  u  is  called  the    'pulse  solution'. 

For  certain  operators  L  it  can  be  proved  that  the  decay  condition  is 
satisfied.     This  is  done  in  the  Appendix, 
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,00 


u  (t,x)e    dt 


t" 


exists  and  is  bounded  independently  of  co  for  oo  >  co  >  C.  Now,  by  integrating 


by  parts  it  is  easy  to  show  that 

T  T 

/_\   iofl"       .      (      icot   ,. 
u(T;e         •  ioj      ue         dt     + 


/ 


ioit  J. 
u.    e         dt  ♦  c 


where  c  is  independent  of  T»     Letting  T  ->oo  we  see  that  the  right  side  of  the 
above  equation  tends  to  a  limit.     Thus  we  may  write 

lim         u(t,x)  e  -     a   . 

t  ->oo 

Put  since  the  integral 


.00 


,.      X   icot    ,. 
u(t,x)e         dt 


converges  it  must  be  that  a  »  0,  We  conclude  that  u(t,x)  ->  0  as  t  ->oo, 
Applying  the  decay  condition,  (117)  may  be  rewritten 

(120)         w(t,x)  «  z(t,x)  ♦  e'^'^  v(x) 


where 


00 

(121)         ^(t,x).  :iil^  -e'^^^  j    u(s,x) 


e    ds 


and 


.00 


(1??)         v(x)  -    u(s,x)  e^"*  ds. 


'o 


It  follows  at  once  that  z(t,x)  -*  0  as  t  ->oo.  The  function  z   is  called  the 
transient  and  v(x)  is  called  the  steady  state  func^tion  of  the  problem  (118), 
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Note  that  v  is  uniquely  determined  by  the  decomposition  (120)  and  the  condition 
that  z  ->  0  as  t  -i>a3.  Furthermore  it  is  easily  seen  that  v  is  independent  of 
the  initial  data  of  the  problem  (118),  for  we  may  modify  that  problem  by  adding 
to  w  the  weak  solution  u'  of  any  problem 

(123)  Lu'  »  0}.   u'(0,x)  -  g'(x)j       (u*  in  T  on  X  -0) 

where  g'  is  of  type  g.  The  resulting  function  w'  «  w+u'  will  be  the  weak 
solution  of  the  problem 

(12U)         LW  -  g(x)e"^'**';    w'(0,x)  -  g'(x)  -  ^  i    (W  in  T  on  x^  -  O). 

But  problem  (12li)  has  a  steady  state  function  which  is  again  given  by  (122), 
The  steady  state  function,  v,  defined  above  is  our  main  object  of 
interest  in  this  paper.  We  shall  attempt  to  show  that  the  formal  expansion  v 
derived  in  Section  1  is  an  asymptotic  expansion  of  v  for  oo  ->oo.  To  this  end 
we  shall  employ  the  integral  expression  (122),  In  Section  6  we  shall  show  that 
V  is  a  solution  of  the  reduced  equation,  Nv  ■  g. 


5,   Asymptotic  expansion  of  the  steady  state  function 

In  Part  I  we  proved  that  under  appropriate  conditions  the  function 
u(t,x)  is  piece-wise  smooth,  i.e.,  u  is  smooth  everywhere  in  KJ    except  on 
certain  hypersurfaces  C°'  where  u  and  its  derivatives  have  jump  discontinuities. 
In  Section  6  we  will  show  that  these  hypersurfaces  can  never  have  a  tangent 
plane  parallel  to  the  t-axis,  so  they  can  be  defined  by  equations  of  the  form 
t  «  T^°'(x),  With  this  information,  the  asymptotic  expansion  of  the  integral 
(122)  is  obtained  simply  by  repeated  integrations  by  parts.  At  each  step  we 
obtain  contributions  from  the  lower  end  point  of  the  integral  and  from  the  dis- 
continuities of  u  and  its  derivatives.  We  shall  assume  that  the  discontinuities 
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of  u  and  its  derivatives  with  respect  to  t  all  occur  on  a  finite  number  of  hyper- 

surfaces  t  -  f   (x)j  a  -  l,.,.,p.  We  write 

^(1)  ^(2)        00 

(125)  v(x)  -  ['  *  [    *  •••  +     u(t,x)e^*dt  . 

Integrating  by  parts 

(126)  v(x)  .  -  ^%S)   -  2:    ?-,_  \_.T-\    \  *■■■*[  <.,(t,x)S^  dt    I    . 


H«re 


0,x)   -  u    (  T^°^-0,x\  . 


(127)  [u]°     .    yx(l^°-^   +0,:^ 
By  repeated  integrations  by  parts  we  obtain 

M    _T  r+1  <,,t(o) 

(128)  v(x)  -  H  (^)    \   u  ^(0,x)  *  YL  e 


T  M+1 


00 


J, 


\(M*1)  ^    ^*   • 


o 

n 


T—      V 

Since  u  is  a  piecewise  smooth  solution  of  Lu  »  u.  +  >   A  u   +  Bu  •=  0,  and 

/?  V»l      V 

u(0,x)  ■  g(x),  we  have  u  (0,x)  =  -  c<  g,  and  in  general 


(129)  u  ^(0,x)  -  (-/)''  g. 

The  operator  oC    is  defined  by  (6), 

Applying  the  decay  condition,  the  integral  in  (128)  may  be  estimated 
independently  of  ca.  We  substitute  an  order  symbol  for  the  last  term  in  the  sum 
together  with  the  integral  term: 
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(130)     .(X) .  -  r  -^  ♦  r  -'"'   '"'  E  (,i)    y  *  °(4i>  • 

r-0  (lo))      a  r"0      l_t_l     « 


If  g(x)  is  sufficiently  smooth,  we  know  from  Part  I  that  u  will  be  smooth  and  the 
terms  involving  [u  ]*^  will  vanish.  Now  compare  (130)  with  (10),  We  see  that 
V  «  y  is  indeed  an  asymptotic  expansion  to  M  terms  of  v.  In  the  general  case, 
where  g  and  its  derivatives  are  discontinuous  across  p ,  we  must  determine  the 
functions  (-1)  )  [u  1'^  which  appear  in  (130),  We  shall  identify  them  with 
the  terms  z^°'^  "^  in  the  formal  expansion  (see  (3I4)  and  (13)).  We  must  also 
identify  the  functions  T^^  appearing  above  with  the  j^'^'   in  (13). 

Thus  we  must  study  the  discontinuities  of  u.  This  is  done  in  the 
next  section. 


6.   The  discontinuities  of  u 

The  present  section  is  essentially  a  repetition  of  Sections  2,  3  and 
6  of  Part  I,  There  is,  however,  one  important  difference;  In  Part  I  we  studied 
discontinuities  in  derivatives  of  u  normal  to  the  discontinuity  surfaces.  Any 
other  non-tangential  direction  of  differentiation  would  do  as  well.  In  parti- 
cular differentiation  in  the  t-direction  is  permissible  provided  this  direction 

it 
is  never  tangential  to  a  discontinuity  surface  ,  In  this  section  we  shall  use 

the  t-direction.  This,  together  with  the  fact  that  the  coefficient  matrices 

A  ,  B  are  independent  of  t, provides  a  considerable  simplification.  Furthermore, 

we  shall  need  the  results  in  terms  of  t-derivatives  in  order  to  determine  the 

#  — 

Discontinuity  surfaces  (i.e.,  characteristic  surfaces)  may  well  have  tangent 

planes  parallel  to  the  t-axis.  This  fact,  which  is  true  for  Maxwell's  equations 
(71)  and  (72),  is  commonly  unrecognized.  However  such  surfaces  will  be  eliminated 
by  imposing  conditions  on  the  initial  data,  g(x). 
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Qu      ]]**  of  (130).     The  calculations  will  not  be  carried  out  fully,  but  only  to 

t^ 
the  point  ^ere  we  may  compare  them  with  Section  1, 

The  function  u  is  defined  by  the  following: 

(131)  I   j  uM  XL  dxdt     -     0  f  cr  all  i\  in    J^   , 

(132)  u(0,x)     -     g(x) 

(133)  u  is  in  -^  , 

and,  for  problems  with  boundaries, 

(13U)         (u  in  T  on  x^  -  0). 

Suppose  u  is  smooth  on  either  side  of  the  smooth  hypersurface  C,  whose 
equation  is  t  •  T(x),  We  introduce  the  regular  transformation 

(135)  ^    '     t  '   Y(x),    y  -  X. 

Here  y  "  y, ,...,y  .  Ve  easily  obtain  the  transformed  operators 

^       V 

(136)  Lu     «     Gu^  ♦  ^     A  u       +  Bu 

^        V"l  ^v 


(137)  MIY    .     -(G/1)^    -f     (A^-TL)       +  B*  XI 

^        v»l  -'v 

where 

(138)  G     -     I  -  ^       T      A*^. 

v-1      *v 

If  we  integrate  by  parts  in  (131)  and  use  the  fact  that  jfL  is  arbitrary  we  obtain 

(139)  gQu]     -    0  . 
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The  symbol  (^u  ]  «  u(Y+0,x)  -  u(T-C,x)  denotes  the  jump  in  u  across  C* 

V  J 

Since  A     and  B  are  independent  of  0,  we  may  differentiate  the  equation  Lu  -  0  to 
obtain  Ln^(^^.i)  "  0  on  both  sides  of  C ,     Hence 

(lUO)  g[u      1+7"     A^fu  /     T  J       +  b|u  /     T  J   «  0:         r  »  1,2,...     . 

Since  lu   I  ^  '^  ^^^   some  r,  we  must  have 
(lUl)         det  G  -  0. 

It  follows  from  (lUl)  that  the  surface  C  is  a  characteristic  surface 

1     k 

for  the  operator  L.  Let  R  ,...,R  be  a  complete  set  of  orthonormal  eigenvectors 

of  G.  We  now  impose  the  multiplicity  condition  of  Section  1.  It  follows  that 

on  C  certain  eigenvalues;  \^,«*.,X  of  G  are  identically  zero  and  all  the  others 

are  non-zero.  Let  R  ,...,r"  denote  the  associated  null  eigenvectors. 

Now  observe  that  u   •  u  ,  and  that  the  function  u   I  is  independent 
f         t^  1-0^  J 

of  0.  Set 

(lUJ)        [u  ^]  -[",]-  (-D'^^'.'^^'Cy)  -  (-l)('*".<'-*^'(x).  r  .  0,1,2., 

z<°'  -  0. 

Then  (139)  and  (lUO)  may  be  written 

(1143)        Gz^^""-^^  -  ofz^^^  '0;  r  -  0,1,2,...  . 

Now  compare  (m3)  with  (l5);  (138)  with  (16);  (lUl)  with  (17).  It  is  clear  that 
we  may  expand  our  z    as  in  (21)  and  obtain  transport  equations  for  the  expansion 

coefficients  which  are  identical  to  the  transport  equations  of  Section  1,  namely 

*  (r) 

(32)  and  (2h)  .  It  follows  that  we  may  identify  our  functions  T  and  z    with 

Some  conftision  might  result  from  our  definition  of  rays.  The  rays  of  Part  II 
are  the  projections  into  the  plane  t  •  0  of  the  rays  of  Part  I, 
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those  of  Section  1  provided  we  can  show  that  they  are  determined  by  the  same 
initial  conditions. 

The  function  u  will  be  discontinuous  across  those  characteristic 
surfaces  C    which  intersect  the  initial  plane  t  ■  0  in  the  manifold  p, 
across  which  g  is  discontinuous.  Let  the  equation  of  C    be  t  ■  T^°'(x). 
Then 

(lUU)         T^"^  -  0   on  r. 

Since  C    must  be  tangent  to  p,  the  projection  into  the  plane  t  ■  0  of  the 
normal  to  C    must  be  parallel  to  the  normal  to  P ,  Thus 

(1U5)         yI°-^     -  p(a)Y^  }        V  -  l,...,n 

V 

Here  Yi»«»»>Y     ^re  components  of  the  unit  normal  to    p   and  p(a)  is  a  propor- 
tionality factor  which  is  determined  by  the  equation 

(lli6)  det  G^°^   -  det(l  -  f*  T^°V)  -  det(T-p(a)E)   -  0,     on   P 


^' 


X 

V 


where 

(1U7)  E    -    r     r,  a", 


v»l 


There  will  be  as  many  characteristic  surfaces  C         emanating  from   p 
as  there  are  distinct  solutions  p(a)   ■  oT^  of  (1U6),  hence  as  many  as  there  are 

distinct  non-zero  eigenvalues  p(a)  of  E.     Precisely  as  in  Section  1,   (page  kh) 

1  V 

we  introduce  the  eigenvectors  B*'  and  J     of  E.     Instead  of  a   •source  condition' 

it  is  more  appropriate  to  define  an 
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Initial  data  condition;   J[g3"0»  ^"  !»•••»?»  °^   P» 
(The  conditions  are,  of  course,  identical  ,) 

By  using  the  fact  that  Lu  »  0,  we  may  obtain  the  values  of  u   on  the 

t 
initial  plane,  t  -  0,  Thus 

(1U8)         u  ^(0,x)  -  (-/)^g(x)}   r  -  0,1,2,...  . 
t'^ 

The  initial  values  for  the  transport  equations  are  obtained  from  the 
fact  that 


(1U9) 


[^  ^(0,x)ir  -  -^e(a)[u  ^J"j  on  T;   r  -  0,1,2, 


. .  •  . 


Here  [  ]'  and  [  ]°'   denote  jumps  across  pand  across  the  characteristic  surface 
C    and  e(a)  ■  t.  1.  The  orientation  of  [  ]*^  is  already  defined.  We  choose  the 
orientation  of  [  ] '  to  coincide  with  the  orientation  of  [  ] '  as  defined  in 
Section  1.  5y  appropriately  indexing  the  characteristic  surfaces  C    to 
correspond  to  the  waves  z    of  Section  1  it  is  not  difficult  to  see  that  the 
numbers  e(a)  will  have  the  same  values  here  as  in  Section  1.  Using  (1U8)  and 
(1U2)  we  may  rewrite  (1U9)  as 

(150)         [/^^~^^g]^^  -E«(=)^^''^^''^   °"  Pi   r-1,2,...  . 

Now  compare  (lUU)  with  (37) J  (lU6)  with  (Ul);  (150)  with  (39).  It  is 
clear  that  the  functions  z^*^^^     and  Y    are  determined  by  the  same  initial 


The  initial  data  function  g  is  identical  with  the  source  function  g  of  Section  1, 
Hence  the  introduction  of  a  'new'  condition  on  g  would  appear  to  be  superfluous. 
However  the  results  of  this  section  are  applicable  to  initial  value  problems  in 
general,  just  as  is  the  whole  of  Part  I.  From  this  point  of  view  we  really  have  a 
condition  on  initial  data.  It  is  this  condition  which  makes  an  analysis  in  terras  of 
t-derivatives  possible.  For  Maxwell's  equations  the  condition  is  particularly  inter- 
esting for,  as  we  have  seen,  it  sheds  new  light  on  the  'divergence  equations'  (73),(7ii), 


-  7U  - 


conditions  as  they  are  in  Section  1.  ffence  the  ftinctions  in  the  two  sections 
are  identical. 

For  problems  with  boundaries  we  must  apply  the  boundary  condition  (13U). 
Suppose  a  characteristic  surface  C   ,  emanating  from  the  manifold  p  ^^  ^^^ 
initial  hyperplane  intersects  the  boundary  B  (x  ■  0)  in  an  (n-1) dimensional 
manifold  ^.  We  now  consider  all  the  characteristic  surfaces  C^  which  inter- 
sect B  in  ^  and  attempt  to  choose  initial  conditions  for  the  transport  equations 
on  these  surfaces  in  such  a  way  that  the  jumps  in  u  and  its  derivatives  will  be 
compatible  with  the  boundary  condition. 

The  analysis  in  Section  1  (beginning  on  page  U7)  is  directly  applicable 
here.  We  need  only  identify  Yx -space  with  tx -space j  'phase  surface'  with  'character- 
istic surface',  etc.  Again  we  require  that  the  shadow  condition  shall  hold.  The 

phase  surfaces  which  are  J ^"coming  I  ^^^^  respect  to  B  must  be  identified  with  the 
^  l^outgoingj 

id  6  ciTQ  366  Si  yv 

characteristic  surfaces  on  which  t  J.        V  as  we  proceed  along  a  ray  from  ^o 

lin creases] 

For  convenience  we  shall  also  apply  the  terras  'incoming'  and  'outgoing'  to  the 
characteristic  surfaces  C^  •  Thus  one  of  the  incoming  characteristic  surfaces 
c'^  is  our  C   .  The  rest  of  them  nay  be  assumed  to  carry  no  discontinuities  • 
The  boundary  condition  may  be  stated  in  the  form 

(151)  b^u  ■  Oj     on  Xj^  ■  0}     V  -  l,...,m. 

The  vectors  b^  are  independent  of  t.     Differentiating  repeatedly  with  respect  to 
t,  we  obtain 

(152)  bu       -O;  onx^.Oi  ^  -  0,1,2;...     . 
Now,  in  obvious  notation, 


If,  by  coincidence,  two  of  the  C         emanating  from   P  intersect  B  in  the  same 
^  ,  we  may  still  treat  them  separately. 
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Using  the  notation  defined  in  (m2),  (153)  way  be  written 

(15U)        i:b^z(p>('>-oi      onBi  ;::j*i--» 

We  observe  that  (15U)  is  identical  with  (63).  It  follows  that  the  z^^'^^^   of 
this  section  and  the  corresponding  functions  of  Section  1  have  the  same  initial 
values  on  B,  hence  are  identical. 

Having  determined  all  the  discontinuity  surfaces  of  u  and  all  the  dis- 
continuities of  u   (r  ■  0,1,2,...)  on  these  surfaces,  we  renumber  them  with  the 

t 
index  a.  Now,  using  the  notation  of  (1U2)  we  may  rewrite  (130)  as 

r-0  (io))      a  r«l  (ioo)         w 

Comparing  this  equation  with  (  70  )  we  see  that  under  appropriate 
conditions,  the  formal  expansion  v  is  indeed  an  asymptotic  expansion  to  M 
terms  of  the  steady  state  function  v.  These  coiditions  will  be  sununairlaed 
in  Section  7. 

We  are  now  in  a  position  to  show  that  the  steady  state  function  v  is 
a  solution  of  the  reduced  equation.  By  this  we  mean  that  for  any  point  x,  not  on 
P,  v  is  in  C  in  a  neighborhood  of  x  and  satisfies  the  equation 


n 
(156)        Nv  -  7~  a\   ♦  (B-ico)v  -  g. 

v«l    *v 


To  show  this,  we  use  the  integral  expression  (122),  and  split  up  the  integral 

at  a  typical  discontinuity  surface  t  ■  T(x) : 

T(x)         00 
(157)  v(x)     -  +1       u(t,x)e^''*dt. 

^o  4(x) 

We  assume  that  they  are  finite  in  number. 
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Differentiating 

(158)  v^       '     \     ^^  e'^'^dt  -e^'^^T^    Vul 

o 

By  integration  by  parts, 

»  00 


(159) 


f     u^e^^^dt     -     -e^^fu]    -u(0,x)   -ico   fue^^dt. 


Using  (158)  and  (159)  we  hav« 


00. 

e      at  +  e 


7"*  A  u       +  Bu  +  u.  , 


1   -JZ     T     A^y  [u1  ♦  u(0,x) 


v«l     *v 


Lu    y  •       dt  +  e  '^     *^(I"J   ■*■  ?(*)• 


But  Lu  -  0  for  t  /  T(x)  and  froni  (139)  G|3"D  "0.     Hence 

(161)  Nv     «     g  for  X  not  on   p. 

Since  g  is  discontinuous  on  p,  it  is  too  much  to  expect  that  v  should 
be  continuously  differentiate  on  p ,  However  it  is  easy  to  see  that  v  is  every- 
where continuous .  This  follows  from  the  integral  expression  (122)  and  the  fact 

■a- 
that  the  discontinuity  surfaces  of  u  are  never  tangential  to  the  t-axis  • 

■5 ^ ■ ■ —- — ^ 

This  last  assertion  is  true  because  we  were  able  to  carry  out  the  construction 

of  the  discontinuities  of  u  without  using  characteristic  surfaces  tangential 

to  the  t-axis.  This  was  possible  only  because  of  the  initial  data  condition. 

Had  we  used  the  formulation  of  Part  I  we  would  have  discovered  that  the  initial 

conditions  for  the  transport  equations  on  such  surfaces  were  zero. 
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7»       SuBimary 

In  this  section  we  suiranarize  the  main  results  of  Part  II.     Statements 
which  refer  to  problems  with  boundaries  are  enclosed  in  parenthesis: 

Giv«n  the  problem 

Lw  -  g(x)e       **;  w(0,x)  -  h(x);  (w  in  T  on  x     -  O) 

where 

1)  L  is  a  symmetric  hyperbolic  operator  with  coefficients  that  are  smooth 
functions  of  x, 

2)  g(x)  and  h(x)  are  of  type  g 

3)  g(x)  satisfies  the  source  condition 

U)   the  multiplicity  condition  is  satisfied 

5)  the  caustic  condition  is  satisfied 

6)  (the  shadow  condition  is  satisfied) 

Then  there  exists  a  weak  solution  w(t,x)  defined  in  the  region  t  >  0  which  is 
unique  in  X  .  If»  in  addition, 

7)  the  decay  condition  is  satisfied 
then  w  has  the  decomposition 

w(t,x)  ■  z(t,x)  +  e     v(x) 

where  z(t,x)  ->  0  as  t  ->oo.  The  steady  state  function  v(x)  is  uniquely  deter- 
mined by  this  decomposition  and  is  independent  of  the  choice  of  the  function  h(x). 
It  is  continuous,  and  except  on  the  discontinuities  of  g(x)  is  a  solution  of  the 
reduced  equation 

Nv  ■  g. 
The  formal  expansion  v  derived  in  Section  1  is  an  asymptotic  expansion  to  M  terms 
(for  o)  ->co)  of  V,  M  will  be  as  large  as  we  please  if  the  coefficients  of  L,  (the 
boxmdary  condition),  and,  except  on  the  manifold  P,  the  function  g(x)  are  suffi- 
ciently smooth. 
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Appendix 

The  Behavior  of  u  for  Large  t 

In  Section  U  we  introduced  the  decay  condition.     This  condition,  un- 
like the  others,   cannot  be  verified  simply  by  examination  of  the  operator  L,  or 
by  a  study  of  the  geometry  of  the  rays  of  the  problem.     Thus,  without  further 
investigation,  we   cannot  even  be  sure  that  the  class  of  operators  which  satisfy 
the  decay  condition  is  not  empty.     Unfcrtmiately,   the  mathematical  literature 
on  the  subject  of  the  behavior  for  large  t  of  solutions  of  hyperbolic  equations 
appears  to  be  practically  non-existent.     In  this  section,  we  shall  attempt  to 
examine  the   subject. 

Our  approach  to  the  problem  is  by  the  well-known  method  of  energy  tn- 
equalities   •     Unfortunately,  we  shall  have  to  place  severe  restrictions  on  the 
coefficient  matrices  of  the  operator  L  in  order  to  arrive  at  the  desired  result. 
It  appears  likely  that  these  restrictions  are  not  inherent  to  the  problem  itself, 
but  merely  reflect  the  inadequacy  of  the  method.     Perhaps  this  discussion  will 

stimulate  further  interest  in  the  question. 

2 
The  energy  method  leads  to  an  estimate  of  the  L     norm  of  the  solution 

u.     Thus,  if  we  introduce  the  function 


/  j  u2(t,: 


(162)  N(t)     -    /      I  u'(t,x)dx       , 

we  can,  under  appropriate  conditions,  show  that  N(t)  ->  0  as  t  -^oo.  Unfortunate- 
ly, this  does  not  suffice  to  guarantee  pointwise  convergence  of  u  to  zero  as 

2 
t  ->CD.  However,  convergence  in  the  L  norm  of  u  and  certain  of  its  derivatives 

is  a  sufficient  condition  for  pointwise  convergence  of  u.  This  is  a  consequence 

of  the  following 

This  method  plays  an  important  role  in  vmiqueness  and  existence  theorems  for 
hyperbolic  equations.  See,  e.g.,  [UJ ,  [$J . 
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Lemma  (Sobolev)  <-  -' :  Let  v(x)  be  a  vector  function  of  x  • 


^ 


» * *•*  n* 


Let  V  and  its  derivatives  up  to  order  n  be  square  integrable  over  the 
entire  x-space.  Then 


(163) 


+  •••  +  v 


fv  dxj  and    jvj   ■  max|v,|,  where  v.    denotes  a  component  of  v. 

The  terms  on  the  right  side  of  the  inequality  include  all  derivatives  of  v  up 
to  order  n,  excluding  those  which  involve  differentiation  with  respect  to  any 
of  the  X.  more  than  once. 


Here    ||v|r 


J 


We  shall  write  the  equation  for  u  in  the  form 


(16M 


Lu  ■  u.    +  A  u       +  Bu     •     0, 


Here  we  employ  the  summation  convention.     In  order  to  apply  Sobolev's  lemma,  we 
must  first  enlarge  the   system  (16U)  to  include  not  only  u,  but  also  the  appro- 
priate derivatives  of  u.     This  is  done,  simply  by  differentiation  of  (I6I4)  with 
respect  to  the  required  x  variables.     By  introducing  vectors  and  matrices  with 
many  more  components  we  may  express  the  enlarged  system  in  vector  form.     Thus 
we  introduce  the  vectors 


»—    u 


(165) 


u   ■ 


n 


u 


V2 


g 


x,x^.»»x 

_      1  t         n. 


^^ 


1 


n 


^1*2 


X   X    . • .x 

U     12         n 
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The  enlarged  systen  new  takes  the  form 


Lu  ■  u.  +A  u   +Bu  -  0, 


It  is  tedious  to  calculate  the  matrix  B  explicitly, 

'"    A^       0      .      .       0         0 


However  it  is  easy  to  see  that 


(166) 


A^ 


0 


0         0 

Lo       0 


A 
0 


Thus  the  k  x  k  submatrices  appearing  along  the  diagonal  of  A  are  all 
identical  to  A  ,  and  the  other  entries  of  A  are  zero.  For  our  purposes  it  is 
sufficient  to  observe  that  the  top  row  of  k  x  k  submatrices  of  the  matrix  B  con- 
sists of  the  matrix  B  on  the  left  followed  by  zero  matrices.  This  too  is  easily 
seen.  We  now  define  a  space  i/     of  testing  function  X\_(t,x),  The  space  v 
is  defined  exactly  as  is  the  space  v     in  Part  I,  except  that  the  vectors  S\. 

Ait 

have  as  many  components  as  does  u.     It  is  convenient  to  introduce  a  special  type 
of  element  of    >f    .     Let 


(167) 


A 


0    J 


Here  Sh  is  any  element  of  M    ,     As  before,  we  introduce  the  adjoint  operator 


M  _fV     -     -  (   a"  JV)  ^     +  B'X^  . 

V 

We  say  that  u  is  a  weak  solution  of  the  initial-value  problem 


(168) 


u 


u(0,x)  .  g(x) 


if  and  only  if 
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(169) 


u  M  X^  dxdt  «  0       for  all  X\_  in  yf  ;       u(0,x)   -  g(x). 


Ve  assurae  that  g(x)  is  of  type  g  and  vanishes  identically  outside  a  compact  region, 
and  that  L  satisfies  the  mviltiplicity  condition.     We  also  assume  that  the  caustic 
condition  is  satisfied  for  the  problem  Lu  ■  Oj  u(0,x)  ■  g(x).     It  is  then  easy 
to  show  that  g(x)  is  of  type  g,  and  vanishes  outside  a  compact  region,  that  L 
satisfies  the  multiplicity  condition,  and  that  the   caustic  condition  is  satisfied 

A/ 

for  the  problem  (168).  It  is  obvious  that  L  is  a  symmetric  hyperbolic  operator. 
It  follows  frcm  Theorem  2  of  Part  I  that  there  exists  a  unique  solution  u  of 
(168)  in  the  class  %    of  piecewise  smooth  functions,  u  is  defined  for  all 
t  >  0. 

We  observe  from  (166),  (167)  and  our  observations  about  B  that 

(17C)         u  M  fv  -  u  M  _n_. 

Here  u  denotes  the  first  k  components  of  u.  Hence  frcm  (169), 

u  M  £]_  dxdt  ■  0   for  all  S\.  in  «  • 

Thus  u  is  the  solution  of  the  problem 

(171)        Lu  =  0;    u(0,x)  -  g(x)j 

>*ich  is  considered  in  the  definition  of  the  dec^  condition. 

Before  deriving  the  energy  inequality,  we  wish  to  recall  the  first 
steps  in  the  calculation  of  the  discontinuities  of  a  solution  of  a  discontinuous 
initial  value  problem   .     By  integration  by  parts. 


(172)  0  -  u  M  _r\.  dxdt  - 


j\  L  u  dxdt  + 


4 


n 


ru    G  u     ds. 


This  is  mentioned  in  Part  II,  Section  6,  but  is  carried  out  in  detail  in 
Part  I,  Section  2. 
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Here 


(173) 


G     -     I    -  T       A 


and    ff  J   "  f(T+0,x)   -  f(T-C,x)   =  T^  -  f_  denotes  the  jump  across  the  character- 
istic  surface,  C,  whose  equation  is  t  «  T(x).     It  follows  that 


(17U) 


Qg  u]  -    gQ^3  -    0. 


From  (17b)  we  observe  that 


(175)      uGu   •  u  G  u^-u^G  u_+  u_G  u^-u_G  u_  •  u^G[_uJ  +  u_G[_uJ  ■  0 


Now,  by  integration  by  parts. 


(176)    0 


u 


L  u  V  dxdt 


T  r  , 

u  M  u  dxdt  -  21   p  G  ujds  +   u 


'^      o  • 


The  C^  denote  the  discontinuity  surfaces  of  u.  We  intrcduce  the  function 
(177) 


N^(t) 


u  (t,x)dx 


and  rewrite  the  adjoint  operator  in  the  form 


(178) 

Here 

(179) 


Ki^  -   -I-_a.-QXi_ 


Q  -  A   -  B  -  B'. 


Now  from  (176),  using  (175),  (177),  (175), 

T 
(180)         0  -  -  I 

o 


a  Q  u  dxdt  +  N^(T)  -  N^(0). 


Here  the  integrals  exist  because  for  each  fixed  t,  u(t,x)  vanishes  outside  a 
compact  region.  This  is  a  conseqaenoe  of  the  same  property  of  g(x). 
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By  differentiation  with  respect  to  T, 
(lei)         2N  N(t)  -   I  u  Q  u  dx. 

We  now  impose  the  following  restriction  on  the  operator  L: 
energy  condition;  For  every  v,  v  Q  v  <  -2y  v  ;  y  >  0;  y  independent  of  x. 
Thus  we  require  that  the  energy  matrix,  Q,  be  negative  definite.  Since  this 
matrix  can  be  calculated,  the  energy  condition  is  an  algebraic  condition  on 
the  coefficient  matrices  of  the  operator  L.  For  a  given  operator  one  can 
easily  determine  whether  or  not  the  energy  condition  is  satisfied. 

Now  (181)  may  be  rewritten 

(182)  N  <  -  Y  N. 

It  follows  that 

(183)  N(t)  <  C^e'^'^. 

Here 

(18U)         C^  -  N^(0)  -   f  S^(0,x)dx  .  (   g^(x)dx. 

From  Sobolev's  lemma, 

(18^)         |u(t,x)  1^  <  2°  f  u2(t,x)dx  -  2"  N^(t)  <  2"cje"^^^  . 

Thus  lu(t,x)  I  decays  exponentially  as  t  ->oo  for  every  x.  This  certainly 
suffices  to  guarantee  that  u  satisfies  the  decay  condition. 
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